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ABSTRACT 

The digital computer simulation of the six degrees of 
freedom equations of motion for the XR-3 captive air bubble 
testcraft is presented. The origin of this computer program 
is the SES Loads and Motion Program developed by Oceanics Inc. 
for the Bell 100 ton (100B) surface effect ship. Modifi- 
cations and procedures used in the revision of the subrou- 
tines to convert the L & M Program from the 100B model to 
the XR-3 model is documented. Measurement data from XR-3 
test runs are used to verify the model for steady-state 
operating conditions in calm water. Computer output for 
turn maneuvers in calm waters and for regular sea conditions 
are included. 



2 



TABLE OF CONTENTS 



X. INTRODUCTION 13 

A. BACKGROUND 13 

B. OBJECTIVES 15 

II. EQUATIONS OF MOTION- 17 

III. PROGRAM SUBROUTINE MODIFICATIONS 30 

A. MAIN 30 

B. SUBROUTINE BOW SEAL 30 

C. SUBROUTINE COLFIL 41 

D. SUBROUTINE INCON 42 

E. SUBROUTINE PROP 50 

F. SUBROUTINE RHS - 56 

G. SUBROUTINE RUDDER 56 

H. SUBROUTINE SIDEWALL 59 

1. Crossflow Drag Terms 60 

2. Slender Body Theory Terms 61 

I. SUBROUTINE STERN SEAL 64 

J. COMMON BLOCK TIE MAP 66 

IV. PROCEEDURES 70 

A. OBTAINING STEADY STATE 70 

B. DIFFERENCES ENCOUNTERED IN THE MODEL 71 

1. Seal Configuration 71 

2. Plenum Dimensions 72 

C. SENSITIVITY TO CHANGES IN INPUT DATA 72 

D. POSSIBLE ERRORS IN THE PROGRAM 73 



3 



V. PROGRAM VERIFICATION 

A. INSTRUMENTATION AND RECORDED DATA 

B. LOADS AND MOTIONS PROGRAM DEVELOPED DATA 

1. Steady State Trim Conditions 

2. Presentation of Unverified Program Data — 

VI. DISCUSSION AND EVALUATION OF SIMULATION RESULTS — 

A. DISCUSSION OF PRODUCT MOMENT OF INERTIA, I - 

B. 5 DEGREE TURNS WITH AND WITHOUT WAVES 

C. TURNS WITH 5, 15, AND 20 DEGREE RUDDER 

D. COMPARISON OF SIMULATION AND MEASURED DATA 

VII. CONCLUSIONS AND RECOMMENDATIONS 

APPENDIX A SES MOTIONS AND LOADS PROGRAM USERS MANUAL — 

APPENDIX B FORTRAN PROGRAM LISTING 

APPENDIX C SAMPLE OUTPUT 

BIBLIOGRAPHY 

INITIAL DISTRIBUTION LIST 

FORM DD 14.73 



76 

76 

79 

79 

82 

95 

95 

95 

99 

99 

109 

111 

183 

234 

239 

241 

243 



4 



LIST OF TABLES 

I. Test Craft Measured Data -- 78 

II. L & M Program Initial Conditions 80 



5 



LIST OF DRAWINGS 



1. S.E.S. Plane of Symmetry 18 

2. Relative Orientation of Reference Frames 22 

3. Static Forces Acting on a Surface Effects Ship 24 

4. Seal Sections Port Half Bow and Stern 31 

5. Fan Flow Characteristics 32 

6 . Ducting Configuration 33 

7. S.E.S. 100 Stern Seal 35 

8 . Bow Seal Hinge Arrangement 37 

9. Bow Seal Configuration 38 

10. Rudder Geometry 47 

11. Propeller Thrust Diagram 53 

12. Resultant Propeller Forces 54 

13. Stern Seal Configuration 65 

14. Thrust Versus Velocity 83 

15. Pitch Angle Versus Velocity 84 

16. Draft Versus Velocity 85 

17. Thrust Change: Speed Vs. Time 86 

18. Thrust Change: Pitch Angle Vs. Time 87 

19. Thrust Change: Z Displacement Vs. Time 8 8 

20. 20 Knot Turn: Rudder Angle Vs. Time 90 

21. 20 Knot Turn: Pitch Angle Vs. Time 91 

22. 20 Knot Turn: Roll Angle Vs. Time 92 

23. 20 Knot Turn: X Y Displacement 93 

24. Five, Fifteen and Twenty Degree Turns: Roll Angle Vs. 96 

Time, I x2 equal to zero. 



6 



25. 



Five, Fifteen and Twenty Degree Turns: Pitch Angle Vs. 97 



Time, I equal to zero. 

26. Five, Fifteen and Twenty Degree Turns: 
Time, I xz equal to zero. 

27. Five Degree Turn: Roll Angle Vs. Time ‘ 

28. Five Degree Turn: Pitch Angle Vs. Time 

waves . 

29. Five Degree Turn: Wave Height Vs. Time 

30. Five Degree Turn: X Displacement Vs. Y 

31. Five, Fifteen and Twenty Degree Turns: 
Time 

32. Five, Fifteen and Twenty Degree Turns: 
Time 

33. Five, Fifteen and Twenty Degree Turns: 
Time 

34. Pitch and Roll Angle: Measured Data 



Yaw Angle Vs. - 98 

100 

with and without 101 

102 

Displacement “ 103 
Roll Angle Vs.’ 104 

Pitch Angle Vs. 105 

Yaw Angle Vs. 105 



107 



7 



TABLE OF SYMBOLS AND ABBREVIATIONS 



AEROD 

BOWSL 

FORIT 

FXAED 

FXBS 

FXPWAV 

FXRUD 

FXSS 

FXSW 

*xx 

P XZ 

I yy 

P ZZ 

INCON 

INTGRL 

K 

'N- 

m 

m b 

N 

P=p 

p b 

PHI=0 

PROP 



aerodynamics 
bow seal 
fourier analysis 

aerodynamic force in x direction 
bow seal force in x direction 

bubble-type force in x direction due to drag 

between bubble and wave surface 

rudder force in x direction 

stern seal force in x direction 

sidewall force in x direction 

mass moment of inertia about x-axis 

mass moment of inertial about xz-axis 

mass moment of inertia about y-axis 

mass moment of inertia about z-axis 

initial conditon 

integral 

roll moment 

length of craft (24 feet) 
mass 

mass of the air bubble 
pitch moment 

rotational velocity about x-axis 
pressure in the bubble 
roll angle 
propeller 



8 



NSRDC 


Naval Ship Research and Development Center 


PSF 


pounds per square foot 


psi=y = 


yaw angle 


Q 


rotational velocity about y-axis 


Qin 


air flow rate in 


Qout 


air flow rate out 


R 


rotational velocity about z-axis 


RHS 


right hand side 


SAM 


shear and moment 


SIDWL 


sidewall 


STNSL 


stern seal 


II 


standard atmospheric reference 


T 


arbitrary steady-state time 


THETA=0 = 


roll angle 


U 


longitudinal velocity 


V 


lateral velocity 


W 


• vertical velocity 


X 


horizontal distance in direction of motion 


y 


horizontal distance perpendicular to direction 
of motion 


z 


vertical distance 


C.G. 


center of gravity 


psf = 


pounds per square foot 


SES 


Surface Effect Ship 


ACV 


Air Cushion Vehicle 


CAB 


Captured Air Bubble 



9 



plenum chamber 



hump speed 



cavity beneath the wet deck of the XR-3, 
formed by the two sidewalls and the bow 
and stern seals 

speed at which testcraft overcomes large 
low velocity drag and begins to operate 
on the cushion of captured pressurized 
air 



10 



ACKNOWLEDGEMENT 



The Authors wish to express their sincere appre- 
ciation to Associate Professor Alex Gerba of the U. S. 
Naval Postgraduate School, for the guidance, assistance, 
and continuous encouragement which he provided during 
the pursuit of this study. 

Also acknowledged is the helpful assistance 
rendered by Mr. Hans W. Doelman, Mathematician, of the 
W. R. Church Computer Center, for his valuable assis- 
tance in the programming work. 

The Authors wish to thank the following W. R. 
Church Computer Center operators ; 

E. V. Donnellan (Supervisor) 

M. Anderson 

J. Shaeffer 

K. Bertwell 

D. Goodwin (Supervisor) 

for their continuous support throughout this project. 
Associate Professor D. M. Layton for his active parti- 
cipation and data collection. He and his technical 
assistant, Mr. Michael O'Dell took great pains to 



11 



provide the authors with sufficient data so that some 
meaningful results could be obtained. 

Last, but not least, the authors wish to extend 
their sincere appreciation to their wives, who labored 
many hours over a typewriter in order that these pages 
might be presented. 



12 



I. 



INTRODUCTION 



A . BACKGROUND 

Conventional monohull displacement type ships are drag 
force limited to speeds of approximately 50 knots. A very 
cursory study of the frictional drag component and its di- 
rect dependence on the square of velocity will bear this out. 

It can be shown that the power plant requirements to over- 
come the drag force at the higher speeds requires such large 
plants that the ship weight increases significantly resulting 
in a speed barrier due to drag forces. The desire for speeds 
a,bove the drag barrier speed has led to some new approaches 
in ship design. The Surface Effects Ship is one of these 
approaches . 

Surface Effects Ship is a broad category which includes 
the ACV, Air Cushion Vehicle, GEM, Ground Effects Machine 
and the CAB, Captured Air Bubble. Navy Magazine of May 1967 
has a very informative article on different ACV and GEM 
concepts. However., we shall confine our discussion to the 
CAB type of craft. 

Captured Air Bubble Concept - The captured air bubble 
craft consists of rigid side walls with flexible soft seals 
between the water surface and the ship at bow and stern. 

This provides a chamber into which air is blown and cap- 
tured. The resulting pressure in the plenum provides a 
vertical force sufficient to allow shallow draft operation. 

The shallow draft produces a lower drag and a higher 
speed for the same thrust than would be had for a corresponding 
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displacement hull. The reason for this effect is that the 
CAB ship drag is composed of two distinctly different type 
of forces. The drag resulting from skin friction on the 
immersed surface is directly proportional to velocity squared. 
The drag resulting from the friction of the bubble at the 
air-water interface is inversely proportional to velocity to 
the 1.566 power. Thus the sum of these two drag forces has 
a point at which an increase in speed actually results in a 
reduction in drag. 

With the development of faster submarines and the spec- 
tacular increases in speed of aircraft the Navy is rightfully 
searching for a method of increasing the surface ship speed. 

A special project office, Surface Effects Ship Project Office 
CSESPO) , has been created to specifically look into this 
problem and to evaluate proposals. 

At the present time two 100-ton Captured Air Bubble 
vessels have been constructed. They are the Aerojet General 
100-A and' the Bell Aerospace Systems 100-B. These two craft 
are prototypes of a projected 2000 ton craft. Prior to con- 
struction of the 100-ton crafts, experimental models were 
constructed in the 3-ton range. One such craft, the XR-3 
built by NSRDC, is presently being used in experiments at 
the Naval Postgraduate School. Additionally, Oceanics In- 
corporated was commissioned to develop a digital computer 
simulation of the CAB ship loads and motions. The program 
delivered to NPGS in October 1972, was a model of the 100-B. 
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B. OBJECTIVES 



Development of the equations of motion for the CAB 
ship were made in terms of the various elements contributing 
separate forces to the different degrees of freedom. The 
objectives of the simulation given in reference 9 are here 
restated . 

1. To provide a computerized form of six degree of 
freedom equations for SES craft that will yield time domain 
outputs of motions for the on-bubble mode of operation. 

2. To determine the motions of specified, SES confi- 
gurations with special references to human habitability, 
platform requirements (for military or commercial appli- 
cations) and maneuvering requirements for the on-bubble 
mode of operation. 

3. To establish a mathematical model for predicting 
maximum wave-induced bending moments, both longitudinal and 
transverse, on SES craft to provide techniques that would 
yield sla'mming responses for on-bubble mode of operation. 

4. To determine the stability derivatives in the ver- 
tical plane from measured data of the vertical plane motions 
of existing SES craft in a random seaway. 

The objective of this thesis is to closely examine and 
attempt to improve the mathematical model. Our method is 
one of revising and changing the 100-B model L&M program so 
that it accurately describes the motion of the XR-3 model 
test craft. Program conversion consists of subroutine modi- 
fications. In certain subroutines it was necessary only to 
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input to the program new data, based upon measured data from 
the test craft. In other subroutines the method of calcu- 
lation, coefficients, or entire equations had to be revised. 
The thesis discusses thoroughly these changes and their cor- 
responding rationale. In addition, a comparison of measured 
test craft data with the L&M program output for steady-state, 
calm water operating condition is used to establish the 
accuracy of the simulation. 

A users manual for this Fortran IV program is in Ap- 
pendix A which includes descriptions of the subroutines and 
the method of inputting data. The program is designed to be 
operated in six degrees of freedom. Since many of the 
crafts parameters are input through the use of data cards 
a thorough reading of the subroutine descriptions and their 
corresponding inputs should be made before attempting to 
exercise this program as it presently exists. 



II. EQUATIONS OF MOTION 



INTRODUCTION 

Motion of any vehicle which is free to move in a three 
dimensional space is made up of six degrees of freedom. 

That is, the vehicle (in this case the XR-3) may have trans- 
lational motion in any one or more of three coordinate di- 
rections and it may also have rotational motion about any 
of these coordinate directions. In order to study these 
motions and the control problems of the XR-3, a mathematical 
model consisting of these six equations is required. These 
equations can then be used as the basis of a computer model 
suitable for simulation studies. In general, these above 
equations will be complex and nonlinear. The initial de- 
velopment and computerization of the six degree of freedom 
motion equations was based on the use of model test data, 
limited XR-3 test data, and theory for estimating the co- 
efficients or constituent terms that appear in the equations. 
The inherent motion characteristics of the XR-3 are determined 
by exercizing this computer model. These results are com- 
pared with available test data and finally used to predict 
other motion characteristics not readily obtainable through 
XR-3 craft testing. 

Early analyses of motions within the vertical and la- 
teral planes were considered separately and revealed that 
linear equations provided an acceptable description of 
motions within the vertical plane (the symmetrical plane for 
the XR-3 as shown in Figure 1.) . However, the linear lateral 
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Lateral Plane Contains 




Figure 1 



SES Plane of Symmetry 



plane equations proved to be inadequate and it was clear 
that nonlinear effects were significant in the description 
of asymmetric motions. 

To produce a resultant functional form for the equations 
that will be the most simple form possible in order to achieve 
an ease of computer representation and solution while still 
including the main factors that affect craft motion, a sim- 
plified nonlinear representation was developed. 

COORDINATE SYSTEM 

The equations of motion describe the craft's movement 
in terms of an inertial reference frame in order to determine 
its position in space. Since the XR-3 is expected to have 
small angles of pitch and roll, the main angle indicating 
orientation of the ship relative to the inertial frame is the 
yaw angle. The range of this angle is large in order to pro- 
perly describe' craft turning motion in the horizontal plane. 

A second set of coordinate axes, used to describe craft 
roll, pitch and heave, is the local level reference frame. 

The horizontal plane of this reference frame remains parallel 
to the horizontal inertial reference plane however its origin 
is placed at the craft's center of gravity above the mean 
free surface and its x-axis is in the craft's heading di- 
rection. The transformation between the local level and the 
inertial reference planes is given by 




cos , sin S\ 0 ij 

-sin 'V , cos S'* , 0 1 

0 , 0 , 1 ^1 
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The craft's heave , pitch, and roll motions are related 
to the local level reference plane through a third set of 
coordinate axes called the Body Frame. The Body Frame, also 
called the ship axis system, has its origin at the ship's 
center of gravity and its positive x-axis is also aligned 
in the direction of the ship's bow, the y-axis is positive 
to starboard, and the z-axis is positive downward. The 
transformation between the body axes and the local level 
axes is given by 
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Finally, to describe the orientation of the fixed body 
frame with respect to the inertial frame a set of single 
axis rotations are performed. These rotations, through what 
are known as Euler Angles, bring the inertial frame into 
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alignment with the body frame. Taking the sequence to be 
positive rotations about the , * and 1 axes, in that order, 
the corresponding rotations are 4^ (yaw) , 0 (pitch), and 0 (roll), 
and the transformation between reference frames is given in 
term S' of these Euler angles by: 
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1 B 
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where : 
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XX 
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COS0COS 4^ 


B 

xy 
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cos0sin V 


B xz 
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sin0 


B 

yx 
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sin0sin0cos 4* - cos^sinM^ 


B 

yy 
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sinj2(sin0sin'V + cos0cos 4^ 


B yz 


•= 


sin0cos0 


• B zx 




cos0sin0cos 4*+ sinJJsin^ 


' B zy 


= 


cos0sin0sin4^ - sinj?cos v { / 


B zz 


= 


cos0cos0 



Figure 2 illustrates the relative orientation of the above 
mentioned reference frames. 

During normal operation of an SES vehicle, roll angle 0 
and pitch angle 0 remain within a range where no great error 
is introduced by using linear approximations for the tri- 
gonometric functions. The resultant transformation between 
the inertial frame and the body frame now becomes: 
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Figure 2 



Relative Orientation of Reference Frames 
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DYNAMICS AND KINEMATICS 



There are three general categories of forces acting on the 
XR-3 type craft. The first category of forces are those acting 
upon the hard structural elements of the ship (rigid body 
forces) such as the sidewalls, rudder, and deck. The second 
group of forces are those acting on the seals and transmitted 
to the ship through their attachments. Finally, the last 
category are those forces which act on the ship and water 
surface due to the air bubble and the lift fan system. Figure 
3 is a simplified illustration of three force categories acting 
on an SES type craft. 

Craft motion occurs as a result of a combination of all 
three force categories. When the vector sum of all these forces 
are in balance, the craft will continue in a state of uniform 
motion. Motions in the plane of symmetry are: forward and aft 

motion or surge; vertical motion or heave; and pitch motion or 
rotation in the vertical plane. Motions in the lateral plane 
are; sway motion or lateral translation; yaw motion a rotation 
about the vertical axis; and roll motion or rotation about the 
forward and aft axis. 

The craft will accelerate in the forward direction as a 
result of propulsive thrust at some equilibrium velocity. The 
pitch motions of the craft are determined by the interactions 
between thrust, drag, lift and gravitational forces. Heave 
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HYDRODYNAMIC 
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c) FORCES DUE TO AIR BUBBLE 



Figure 3. Static Forces Acting on a Surface Effect Ship 
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motions result from variations in bubble pressure, buoyant 
forces on seals, sidewalls, and waves reacting against the 
weight of the craft. 

Motions in the lateral plane result from sideforces 
acting through the center of gravity to produce lateral 
translation or sway, sideforce offset longetudinally from 
the center of gravity to produce angular yawing motion, and 
lift force offset laterally as well as sideforce offset 
vertically from the center of gravity to produce angular 
rolling motion. 

In order to describe the above motions in a simple and 
accurate manner which is conducive to computer programming, 
a correct choice of kinematic variables and coordinate sys- 
tems must be made. When the equations of motion have been 
solved for that set of variables, all other variables can 
then be evaluated using the relationships between kinematic 
variables; and when the orientation of the body frame has 
been established, all other vectors however expressed can be 
resolved into any common frame using the previously given 
transformations between coordinate systems. 

The equations of motion have been based on the following 
chosen kinematic variables: the Euler angle, and their time 

derivatives : vV,e,0; M/,e,0; w, e ,0; and the components of 

the velocity of the origin of the body frame, resolved in the 

direction of the local level frame axes; and the time deri- 
vatives of these velocity components, u L , v-^, w L , and 

• • • ^ 

Vjj, w L , where u L is the component in the i^ direction, 
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v is in the j direction, and w T is in the k_ direction, 

1j ±J J-J J_i 

The equations of motion are to be solved for 

• • t * * * • • • 

V, e, 0 , u L , v L , w L . 



Two advantages to this coice of variables and reference 
frames are; first, the simplicity in the integrations re- 
quired to determine the orientation of the body frame and 
the location of its origin, and second, is that it. facili- 
tates a reduction of degrees of freedom. That is, the e- 
quation for Z can be easily reduced to obtain a representation 
of the heave-only mode which is a readily interpretable 
motion . 

One inescapable disadvantage to the above choice is 
that the transformations required to resolve components into 
local level frame axes and to express relationships in terms 
of selected kinematic variables are nonlinear. In most 
cases this process has been simplified by using only the 
first order terms in the Euler angles. It is to be noted 
that, in comparison to the bubble lift force, the other 
forces acting on the craft are generally small; hence, first 
order transformations from body frame to local level frame 
axes actually provide representative second order accuracy. 

As previously stated a net unbalance of forces or mo- 
ments along or about any axis will cause the craft to ac- 
celerate along or about that axis. This acceleration is 
called a kinetic reaction. This principle known as D'Alemdert's 
principle can be more formally stated that when the kinetic 
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reaction contributions are included, the sum of all forces 
and moments vanish. Stated mathematically we have: 

0 = X -^F. i 

1 l 

o = y -2;f. * 

l l J 

A 

0 = z -?f. k 

1 1 

0> = K - 2.M . l 

1 1 

0 = M - ?-M . t' 

i i J 

0. = N -Zm. k 

1 1 



where X,Y and Z are the kinetic reactions and F^ i“, F^ '5 

A A A A 

and F. k and M. 1 , M. i , and M. k are the various forces 
i i i i 

and moments due to thrust, drag, gravity, bubble pressure, 
and buoyancy in the X,Y and Z directions of the local level 
reference frames. 

To evaluate the kinetic reactions it is assumed that 
the moments of inertia of the craft are constant with re- 
spect to the local level coordinate system. Now, the direct 
application of Newton's laws of motion for a rigid body re- 
lative to this coordinate system leads to the following 
equations of motion or kinetic reactions: 

m (u - v*V) = X 

m (v + u'V) = Y 

m w = Z 



A. 



*xx 0 - W- ^y^ 9 " K 



Jyy 9 - J x 2 W + I- 
t z 2 r - *x 2 f> - N 



XX 



= M 



27 



/ 



Computations using analog computers have revealed that the 
nonlinear inertial coupling terms included in the above e- 
quations contributed only negligible differences to the re- 
sponses of expected forcing functions (see reference 9) . 

Hence, the decision was made to simplify the equations of 
motions to be used by deleting all nonlinear inertial coupling 
terms. The resulting equations are then given by: 



mu = X 

mv = Y - mu*V 

mw = Z 



• * 

I 0 

xx 



-I 

xz ' 



i e 
yy 

tzz*- 



= M 



I 0 

xz 



K 



N 



Note: The u *4‘ / product is retained since it is not a small 

quantity. 

The above simplified equations can be represented in 
matrix form by: 
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xz 
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N 



This matrix representation is used to establish a procedure 
by which the digital computer solves the equations of motion. 
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Detailed descriptions of the component forces and mo- 
ments which contribute to the motion of the XR-3 craft are 
provided in the subroutine descriptions that follow this 
section. 
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III. PROGRAM SUBROUTINE MODIFICATIONS 



A. MAIN 

The only changes to the main portion of the program was 
the introduction of the term VEL and the corresponding for- 
mat change so that velocity is written in knots instead of 
feet per second. 

B. SUBROUTINE BOW SEAL 

The XR-3 is presently equipped with soft air-spring type 
bow and stern seals. These seals were designed and cons- 
tructed by NSRDC and installed by NPGS . The seals illustrated 
in Figure 4 are double-celled air bags consisting of a 46" X 
120" frame and fabric ends, sides, and bottom. Twelve e- 
qually spaced steel spring stiffeners reinforce the bottom 
of the seals.' The divider membrane has holes in it which 
serve as a damping device. Winches and cables have been 
attached to the rings shewn in Figure 4 , so that the shape 
of the inflated seal may be changed. 

To inflate the seals ducted fans are used. The char- 
acteristics of these fans are shown in Figure 5 and the 
ducting arrangement in Figure 6. As can be seen from 
Figures 5 and 6, two ducted fans supply air to the bow seal. 
The rate of flow is 125 cubic feet per second per fan; hence, 
the total flow to the bow seal is 250 cubic feet per second. 
The return ducting, shown in Figure 6, allows the vane axial 
fans to deliver flow and not cause the seal to inflate like 
a balloon. There is sufficient fan flow to produce a pressure 
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Figure 
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FAN FLOW CHARACTERISTICS 
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DUCTING CONFIGURATION 



differential between the seals and the plenum chamber of 
about 1 to 2 psf. depending on the plenum pressure. 

The stern seal of the 100B is of similar construction 
to the bow seal presently installed in the XR-3 and there- 
fore the 100B Stern Seal Subroutine was used as the basis 
for computer representation. By comparison of Figures 4 and 
7, the XR-3 seal configuration and the 100B stern seal con- 
figuration respectively, we can note strong similarities. 

Both are of pressurized fabric-type material that is assumed 
to have no inertia or dynamics and are constructed with a 
vented membrane separating the lobed air bags. The throttle 
shown in Figure 7 is merely a controllable path for the air 
to return to the plenum. This throttle corresponds to the 
return ducting previously mentioned. 

The process of modification of the program involved 
the identification of the elements of the 100B stern seal 
calculations so that proper subsititutions could be made 
without destroying the continuity as well as the way in which 
the program coupled with itself. 

The first computation considered was the calculation of 
the bow seal gap. GAP is here defined as that space resul- 
ting from lifting the seal beyond the surface of the water, 
that is, an air gap resulting in a loss of air from the 
plenum and increased air flow. The equation is here repro- 
duced without its iterative subscripts. 
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SES TOO STERN SEAL 




Figure . 7 



GAP = - (Z+ZS-XX (3 , K) * THETA +YY ( 3 , K) *PHI+ETA (3 , K) ) 

XX (3, K) = XBS+RBS*COS (ANG) 

YY ( 3 , K) = RBS*SIN (ANG) 

This calculation results in the distance DS as shown in 
Figure 9. 

The purpose of the XX(3,K) and YY(3,K) calculations are 
to locate points along the perimeter of the hinge. The dis- 
tances XBS, XX(3,K) and YY(3,K) can all be found in Sub- 
routine INCON. RBS is not required for the XR-3 seal and 
is removed. The reason for the calculation of RBS was that 
the 100B hinge is semicircular in shape as shown in Figure 8. 
The XR-3, on the other hand, has a seal that has its hinge 
on the bow of the ship and is affixed straight across, ath- 
wartships, therefore XX (3, J) = XBS. The distance, YY(3,J) 
is an incremental distance athwartships and is also cal- 
culated in INCON. 

• YY ( 3 , J) = 0 . 5*XBBW+ ( J-l) *DELYBS 
In this way we have located a series of straight line seg- 
ments running athwartships at the seal hinge location. 

Initially the overhang of the bag portion of the 100B 
bowseal had been calculated as ELBAG by Subroutine FGl. 

Since the XR-3 bow seal has no overhang, this calculation 
has been removed. The bow and stern seal of the XR-3 are 
basically alike and therefore the only differences are in 
the definition of terms used in the equations. The cor- 
responding changes have been made in INCON, RUDDER, RHS , 
and wherever else the Common Block ties were made. New 
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BOW SEAL CONFIGURATION 



terms that were included match their corresponding terms in 
the 100B Stern Seal Subroutine. With reference to Figure 9, 
theseo are XBF which is the length of the seal from hinge 
to transom, THBSI which is the angle between the craft ver- 
tical and the seal, BLEAK which is the seal base leakage area, 
ZBS which is the distance from the CG to the seal hinge, 

COSBS and SINBS, which correspond to SINTH and COSTH in the 
100B program. These later two terms are the SIN and COS of 
the THBSI bow seal. 

The following is an equation by equation discussion of 
the XR-3 Bow Seal Subroutine. The first few statements serve 
to zero the force, GAP vector, SKI vector and leakage area. 

(1) DELPBG = PBS - PB 

Equation (1) calculates the differential in pressure 
between the bubble (plenum) PB, and the bow seal bags, PBS. 

PB is calcualt'ed in Subroutine RHS and PBS is an INCON Sub- 
routine input. 

(2) FBAR = PB - PINF 

(3) DELP = PBAR 

(4) IF (DELP. LT. 0.0.) DELP =0.0 

Equations (2) through (4) calculate the pressure dif- 
ference between the bubble, PB, and the sea level value PINF, 
insuring that DELP is never negative since we cannot create 
a vacuum in the plenum. 

(5) SINDIF = SINBS - COSBS * THETA 

(6) COSDIF = COSBS - SINBS * THETA 

COSBS and SINBS are calculated in Subroutine INCON. 
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THETA is the pitch angle and is considered small enough to 
use the small angle approximation for the sin. Therefore, 
equations (5) and (6) represent adjustments to the INCON 
values, COSBS and SINBS, due to pitch angle. 

(7) XI = XBS+ZBS* THETA -XBF*SINDIF 

(8) Z1 = -Z-ZBS+XBS*THETA-XBF*COSDIF 

Equations (7) and (8) locate a line on the water's sur- 
face of the seal based upon the draft and pitch angle. Re- 
ferencing Figure 9 which assumes no pitch or roll it can be 
seen that the distance determined by equation (7) is indeed 
XI. Similarly we see that equation (8) locates the distance 
Zl. It is assumed that the forces acting upon the seal and 
the craft act along this line which is centered in the wetted 
area at the water's surface. The assumption is that the seal, 
being a.‘ buoyant body does not penetrate sufficiently the 
water's surface, but floats upon it. 

The next task is to calculate the wetted area of the 
seal and the gap, if any. The seal may be divided into as 
many as eleven segments according to NSTA(3). 

(9) ELSKI (K) = (ETA ( 3 , K) -DETABX (K) ) * (XX (3,K) -XI) -Z1+ 

W ( 3 , K) *PHI 

(10) GAP (K) = ELSKI (K) 

(11) IF (GAP (K) .LT .0.0) GAP(K)^0.0 

ELSKI represents the wetted length of the small segment 
corrected for waves, (ETA (3 ,K) -DETABX (K) ) : and roll, PHI. 

Equation (11) provides that as long as ELSKI is negative no 
gap exists. Equation (12) computes an average of ELSKI over 
the segments previously determined. 
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(12) ELSKIA = (ELSKICJ+1)+ELSKIA(J) )/2.0 

Obviously if the average wetted length computes to be 
less than zero there is no need to make any of the following 
computations and gap is the only meaningful term. However, 
as is usually the case, ELSKIA has a value and it is then 
necessary to next test ELSKIA and determine if it exceeds 
the maximum wetted length permissable. If so, it is then 
set to the maximum value. This is reasonable since the 
maximum wetted area would exist when the top of the plenum 
was contacted by the water and the computations would stop. 
The succeeding equations calculate a coefficient of drag 
based upon the craft velocity and then the resulting drag 
forces. Similarly the remaining forces and moments are 
calculated . 

The leakage flow, QL, was not permitted in the 100B 
unless the seal broached the surface and GAP became non-zero. 
In the case of the XR-3 it was found that there existed a 
leakage area between the sidewalls and the bow seal itself. 
This gap was calculated and input as CFBS in Block 6 of 
Subroutine INCON. This small area accounts for a flow of 
approximately 10 cubic feet per second. 

Inputs to Subroutine BOWSL will be discussed in the 
section covering Subroutine INCON. 

C. SUBROUTINE COLFIL 

The Colfil Subroutine controls output, both form and 
content. Since substantial changes to the Oceanics L&M 
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program have been made, a user's manual for this subroutine 
is included in Appendix A. The modifications herein des- 
cribed were made for the purpose of arranging the output 
in a more convenient form for the purposes of this report. 
Presently time histories of 16 variables may be printed in 
two tabled summaries of values and up to 10 graphs plotted 
in each run. Additional variables may be output by a single 
change to the present format. The graphs may be drawn on 
either the off-line (calcomp) or on-line (print) plot. 

Changing this routine necessitated minor changes in both the 
RHS and INCON Subroutines. 

The changes in RHS consisted of modifying the method 
of writing the variables on the temporary storage device. 

Only one temporary storage device is used. This was done 
in order that the list can be searched by variable number 
instead of searching by variable name. 

The .changes in INCON necessitated the creation of new 
blocks 20', 21, and 22. Information input and the form of 
that information can be found in the user's manual, Appendix 
A. 

D. SUBROUTINE INCON 

The INCON Subroutine is the most difficult of all the 
subroutines to discuss since it reads and computes the ini- 
tial conditions of all variables and constants which affects 
how the entire program functions. The thorough understanding 
of this subroutines! function and operation are a prerequisite 
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to exercising the program. To assist in this understanding 
a description of the operation of the Subroutine and where 
necessary full description of the important functional com- 
putations is presented in this section. 

The INCON Subroutine is divided into numbered blocks 

which are controlled by the read statement number 10. (see 

Appendix B) . Notice that the method is to read a 3 digit 

*. 

integer followed by a 2 digit integer. The first 3 digits 
are called the block number and the next 2 digits the option 
tag. By executing the computed go to statement according to 
the block number the statement next executed can be controlled. 
Consider as an example the number 01601. The block number 
is 016 and the option tag is 01. Hence, the block number 
will cause the program to flow to statement number 1600 by 
the GO TO control and then be routed., to statement 16 05 by 
option tag number in another GO TO control. Using this flow 
control all of the necessary information is input to the 
program from the card reader and disk. 

What follows is not a bock by block description of the 
input, but rather a description of those places in which 
changes have been made and in which particular calculations 
of interest occur. A block by block explanation of the input 
is contained in the user's manual Appendix A. 

Block 2 option 1 

This input block is the summary mass and inertia input. 
Option 2 allows a discrete mass input to be read. It is 
important to point out the source of the numbers presently 



43 



being used for the XR-3. As a first approximation the moments 

of inertia found on page 1 of reference 2 were scaled up to 

the present XR-3 configuration. The craft weight given in 

reference 2 was 4192 pounds and a ratio, of 5813 = 1.39 was 

4192 

used to compute the following inertias. 

IX = 2941 slug ft 2 
IY = 9320 slug ft 2 
IZ = 11300 slug ft 2 

However, the numbers presently employed are those pro- 
vided by NSRDC and are: 

IX = 2870 slug ft 2 
IY = 9320 slug ft 2 
IZ = 10580 slug ft 2 

As can be seen, the difference is minimal. No information 
was available for the product inertia 1^^ and therefore it was 
assumed to be zero value. 

Block 4 option 2 

The appendages section of the program has been removed 
since there are no fins or other appendages on the XR-3. 

Block 5 

The measured data inputs to the stern seal block con- 
sists of the X and Z coordinate of the seal hinge, the seal 
angle, the length of the seal leading edge, the pressure 
differential between the seal bag and the plenum chamber, 
seal leakage orifice coeffficient and the base leakage area. 

It was assumed that the leakage orifice of the XR-3 had 
the same characteristics as the stex'n seal orifice of the 
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100B. Thus the seal leakage orifice coefficient remains 
unchanged. 

The base leakage area refers to that area over which 
the leakage coefficient exerts influence. This is the gap 
that allows a flow from the plenum to the atmosphere. The 
leakage flow is the control over the plenum pressure, fan 
flow rate and thence the fan horsepower* calculations . The 
leakage area is distributed throughout the craft and it is 
necessary that this number be adjusted to* bring the above 
three calculations into balance at a steady state operating 
condition. If the base leakage area is reduced it causes 
the craft to either rise, or the pressure to increase, or 
the fan flow rate to be reduced. The number presently being 
input is the result of such an adjustment and the three cal- 
culations are in agreement with experimental data concerning 
draft and plenum pressure. 

Block 6 

These bow seal inputs correspond to those of the stern 
seal in Block 5. 

Block 7 option 2 

The critical Froude Number was obtained as follows. 

. . V 13.4 

Fn critical =^/gl"= ^20^ (32) = * 556 

Where: V = speed in feet per second (critical or hump 

speed) 

L = length of the plenum (nominal) in feet 
G = gravitation constant 
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Block. 8 option 1 



Those inputs that correspond to thrust have been re- 
moved and are found in block 16 . 

Block 9 option 1 

Although there is no rudder on the XR-3, it was con- 
sidered necessary to provide information for the rudder 
effect produced by the propulsion system. Refer to Figure 10 
which shows a crossectional view of one of the outboard motors. 
Dimensions in the figure are inches and must be converted to 
feet and square feet for input. 

The first assumption made was that the effect of the 
outboard motors could be separated into two parts, first, 
the forces contributed only by rudder action and secondly the 
forces which result from the thrust. 

The second assumption was that the dashed line represents 
the water surface under most operating conditions. That is, 
at the stern, the water level doesn't change much; at least, 
not enough’ to reflect great rudder force changes. Addi- 
tionally, waves do not interrupt flow at the rudder. There- 
fore, all that is required is to determine the respective 
areas and lengths, convert them to feet and enter them pro- 
perly through the card reader. The wetted surface rudder 
effect area, therefore, is that region which is cross hatched 
in Figure 10. 

Rudder span is defined as the vertical length of the 
rudder. This distance is the wetted surface distance. Based 
upon the foregoing assumptions only 2/3 of the block marked 
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RUDDER GEOMETRY 




Figure 10 



A in Figure 10 through which the water line passes con- 
tributes to wetted area. 

2 

Therefore, rudder span = “(4%)+ll?5 = 14.5" = .(1-21 ft.). 

3 

The rudder aspect ratio is defined as the rudder span 
squared divided by the area, a dimensionless term. 

(1. 21) 2 = 2.15 

.68 

Total rudder area is the summation of the cross-hatched 

area; A and B. 

2 _ 

3A+B = R 

A 

The area of B is taken to be the average width 5.75 
inces multiplied by the length. The area of Block A is 
taken to be a rectangle 11.75 inches by 4.5 inches with a 
small triangle removed. It is felt that these approxi- 
mations will not lead to any gross force inaccuracies. 

= 31.9+66.1 = 98.0 inches^ (.68 feet^) 

Block 10 , 

The aerodynamic inputs are the reference length and 
referenc width. These dimensions are used to compute a 
reference area to be used in conjunction with the parti- 
cular values of the coefficients obtained from the wind 
tunnel tests. 

Block 15 

The shears and moments calculations have been removed 
and replaced by a dummy subroutine SAM, since no stress 
analysis will be performed on the XR-3. This removal 
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provides a convenience for the purposes of this report. 

Block 16 option tag 1 

Entries input in block 16 option tag 1 permit the user 
to input a predetermined time history of thrust on the star- 
board screw as a map. That is, either constant or variable 
thrust may be selected. The user must input data points 
from the map of the variable thruster and the function sub- 
program FG1 is used to determine the functional value of 
thrust at any given time. A maximum of 24 data points may 
be used to describe the map. 

Block 16 option tag 2 

Exactly the same as option tag 1, except that these 
entries apply to the port screw. 

Block 16 option tag 3 

The entries in block 16 option tag 3 allow the user 
to input a rudder motion map. Positive rudder angles are 
right rudder angles which cause the craft to turn to the 
right . 

Block 17 

Option tags_l and 2 permit the user to vary bow and 
stern seal differential pressures as the program run pro- 
gresses. The differential pressure is taken to be that 
between the seal and the plenum. This represents a revision 
to the original program. 

Block 19 option tags 1,2, and 3 

These inputs allow the user to input fan map data and 
thus control air flow in both the seals and plenum. 
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E. SUBROUTINE PROP 



The forces and moments associated with the propulsion 
and steering systems for CAB craft are somewhat unique for 
each propeller and rudder system type. In the original 
Oceanics' program, written to simulate the Bell 100 ton 
SES , the system described consisted of two fixed super- 
cavitating propellers each with a rudder mounted directly 
astern. Both the propellers and rudders were located close 
to the transom of the craft. In the XR-3 , the propulsion 
system consists of two forty horsepower outboard motors each 
is fitted with a two bladed screw. The port screw has a 
right hand rotation and the starboard has a left hand rotation. 
With this system steering is accomplished by rotating these 
motors in the horizontal plane. In turn, the main thrust 
vector produced by each propeller is rotated in the hori- 
zontal plane to provide a desired turning moment. Additional 
turning moment is provided by the motor/propeller housings 
which are essentially shaped like a rudder foil. However 
the surface areas of these housings are significantly re- 
duced from typical rudder areas associated with a fixed 
propeller system. The rudder characteristics are further 
described in the RUDDER Subroutine. 

The rotation of the propeller thrust vector in the 
horizontal plane adds a complication factor to the XR-3 
description which was not present in the 100B craft simu- 
lation. Hence, the original PROP Subroutine required major 
changes in order to make it properly describe the XR-3 craft. 
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In the XR-3 description, as in the 100B, specific pro- 
peller characteristics such as; RPM, Blade pitch, and the 
various propeller coefficients are not included to determine 
the thrust and torque produced by each propeller. In both 
simulations the thrusts are either input values or cal- 
culated from the equations of motion for the associated 
input velocities. 

The propeller representation allows either propeller 
to provide any desired thrust and side thrust independent 
of the opposite propeller. This feature is available 
through the use of function routine FG1. FG1 is used to 
describe the thrust of each propeller as a function of time. 
When the program operates in this mode the velocities and 
accelerations of the XR-3 are obtained from the equations 
of motion as a result of the input propeller thrusts. An- 
other operating mode provides for velocity to be input and 
the resultant total propeller thrust is obtained from the 
motion equations. In this mode of operation, it is assumed 
that the total thrust is divided equally between both pro- 
pellers, and function FG1 is not used since the propeller 
thrusts are not input. 

The above two modes of operation allow for maximum 
program flexibility and simplicity. Test run data of both 
thrust and velocity is available from the installed XR-3 
instrumentation. Thus, the simulation can easily duplicate 
any given set of operating conditions. 
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The force and moment equations for the propeller 
system of the XR-3 are derived from the propeller thrust 
diagram given in Figure 11. The figure displays the main 
thrust and side thrust vectors for both propellers positioned 
with right rudder angle a, positive pitch angle 9, and ne- 
gative roll angle 0 present. The thrust vectors act through 
the port and starboard propeller positions; xp, yp , and zp, 
from the craft center of gravity in the Body Reference co- 
ordinate system. These thrusts are then converted to forces 
in the local level reference coordinate system for the 
calculation of craft pitch, roll, and heave and forces in 
the Inertial Reference frame for calculation of craft tra- 
jectory. That is, the craft yaw angle^, and x and y dis- 
placements . 

The forces associated with the propeller arrangement 
of Figure 11 are depicted vectorially in Figure 12 in the 
Local Level Coordinate system. Employing the simpli- 
fications' used in the derivation of the equations of motion 
namely that the pitch angle 9 and roll angle 9 are limited 
to small values, such that sin0-0, cos0-l, and sin0-0, cos0~l 
the following force equations result: 

T12 = THRUST 1 * COSa 
T2X = THRUST 2 * COSa 
ST1X = STHRUST 1 * SINa 
ST2X = STHRUST 2 * SINa 
TlY = THRUST 1 * SINa 
T2Y = THRUST 2 * SINa 



52 



lAiw \ n icnMM i 




53 



THRUST 2 



RESULTANT PROPELLER FORCES 
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ST1Y 


= S THRUST 1 * COSa 


ST2Y 


= STHRUST 2 * COSa 


T1Z 


= THRUST 1 * COSa * 9 


T2Z 


= THRUST 2 * COSa * 9 


ST1Z 


= STHRUST 1 * COSa * 0 


ST2Z 


= STHRUST 2 * COSa * 0 



By summation of forces in the three coordinate direc- 
tions the following port and starboard propeller and total 
force equations result: 



FXS = 


= T1X + STIX 


FXP = 


= T2Z - ST2X 


FX = 


= FXS + FXP 


FYS = 


= T1Y - STlY 


FYP = 


= T2Y + ST2Y 


FY ' = 


= FYS + FYP 


. FZS = 


= -T1Z - STlZ 


’FZP = 


= -T2Z + ST2Z 


FZ = 


= FZS + FZP 



The moment equations associated with the preceeding force 
equations are given as: 



FKS = 


= FZS * YP - FYS * ZP 


FKP = 


= FZP * YP - FYP *: ZP 


FK = 


= FKS + FKP 


FMS = 


= -FZS * XP + FXS * ZP 


FMP = 


= FZP * XP + FXP * ZP 
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FM = FMS + FMP 



FNS = -FXS * YP + FYS * YP 
FNP = FXP * YP + FYP * XP 
FN = FNS + FNP 

Throughout the preceeding representation it is assumed 
that the magnitudes of the main and side thrust vectors 
remain independent of the various rudder deflection angles a. 
It is also assumed that both port and starboard propellers 
rotate together through the same rudder defelction angles. 

F. SUBROUTINE RHS 

Subroutine RHS contains the expressions of the right 
hand side of the force and moment differential equations. 
These forces are calculated and summed in this subroutine. 

Changes in this subroutine consist of modification of 
the Bubble Drag Coefficient of Friction, Common Block ties, 
and Centep of Pressure Location. The rationale for these 
changes iS presented in the section concerning the presen- 
tation of data. 

G. SUBROUTINE RUDDER 

The rudder subroutine used in the XR-3 simulation re- 
mains essentially unchanged with respect to its force and 
moment equations. However, some changes are introduced in 
the calculation of rudder deflection angle 'a' and in the 
input data to reflect the proper XR-3 rudder location and 
dimensions . 
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Rudder deflection angle can either be input as a con- 
stant angle in the input data set or it can be made to vary 
with time by using function FG1 . Note, this rudder de- 
flection angle is the same angle input to PROP subroutine 
to control the main thrust vector orientation. 

As previously noted in the description of subroutine 
PROP, the motor/propeller housings constitute the XR-3 
rudders. Physical measurements of these housings were used 
as rudder dimensions in the input data set. Each rudder/ 
housing is located directly astern of its associated side- 
wall. They are vertically mounted and rotate in the hori- 
zontal plane. 

The analysis of the XR-3 rudders assumes that they pro- 
duce only lateral and drag forces. Lift or Z direction 
forces are considered negligible and hence they are set 
equal to zero. 

The -development of the force equations associated with 
the rudders of the XR-3 follows the analyses from references 
9 and 7. The lateral force produced by both rudders is 
given by : 

FY=2*RH0*U*U*A r * (1+d/h) *RCLB* (a- (1+d/h) *V h /U) 

where is the rudder area, RHO is the water mass density, 

U is the velocity in the X direction, d is the draft at the 
rudder, h is the draft plus the rudder span or length, RCLB 
is the lift coefficient for the rudder given as 
RCLB = 2AR/ (AR+3) where AR is the aspect ratio, a is the 
rudder deflection angle, is the first order approximation 
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rudder velocity in the lateral direction is given by: 

V, =V+XR*R-ZR*P 
h 

with XR and ZR the coordinate locations of the rudder center 
with respect to the craft center of gravity. 

The drag force produced by the rudders is given by 

FX = -2*CD*Ajj- + RHO*U 2 . 

The drag coefficient is given by 

CD = 2*CFR +PI/8 * (Tr/C) 2 * (1+G (h-d) /U 2 ) + 

RCLB*V//U 

where G is the gravity, Tr is the rudder thickness, C is 
the rudder chord, and CFR is the coefficient of friction 
given by the Schultz-Griinow formula 

CFR = .427/ (Log REY - . 407) 2 ' 64 

where REY • = U*Ar/ (h-d) /ENU . 

REY is commonly referred to as the Reynold's Number and ENU 

♦ * — S 2 

is the kinetic viscosity = 1.28 x 10 ft /sec. 

The two preceeding force equations produce the following 

moments : 

FK = -ZR*FY 
FM = ZR*FX 
FH = XR*FY 
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H. SUBROUTINE SIDEWALL 

Changes made to the sidewall subroutine consist of the 
removal of force computations associated with ventral fins 
which are not present in the XR-3 testcraft. The remainder 
of the subroutine is unchanged from the Bell 100 ton craft 
version described in reference 9 . 

The hydrodynamic and hydrostatic forces and moments 
acting on the sidewalls are found by application of slender 
body theory. This theory, based on motion of a foil in a 
viscous medium, developes methods for finding the lift and 
drag forces associated with moving underwater bodies. These 
methods are primarily linear in nature as shown in refer- 
ences 1 and 9 with the effect of instantaneous changes in im- 
mersed sectional area, draft, etc. due to orientation changes 
arising from heave, pitch, and roll motions reflecting the 
most significant nonlinear variations. The sidewall forces 
also contain certain nonlinear terms due to cross-flow drag, 
which are 'important for the case of very low aspect ratio 
lifting surfaces as are the XR-3 sidewalls. The sidewall 
forces will then be composed of two separate contributions; 
first the terms due to sidewall buoyancy and slender body 
hydrodynamic reactions and second the effect of cross-flow 
drag, i.e. 

R = FH + FD 

sidewalls sw sw 

where FH forces result from slender body theory, and FD forces 
are the result of the cross-flow drag terms. 
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Although no changes were made to the sidewall sub- 
routine force and moment equations they were considered 
of adequate importance to be hete' discussed. 

Before computing the forces on the sidewalls, the 
subroutine first calculates the draft of or the gap be- 
neath the sidewalls. The draft/gap is found for each X- 
direction section after it has been corrected for craft 
roll and pitch angles and waves. If gaps are present, then 
the total sidewall leakage area ALSW is computed by summing 
the products of gap heights and sidewall section lengths 
for both sidewalls. 

1 . Cross-Flow Drag Terms 

The cross-flow drag force components are assumed 
to only be significant for the lateral force, yaw moment, 
and roll moment. While such forces are also present in the 
vertical force and pitch moment, they are not as large re- 
lative to the predominant forces due to bubble pressure and 
are neglected. The lateral force due to cross-flow drag on 
a single sidewall is given by: 

FYD=-RH0/2*CDC*VREL*ABS (VREL) *DELX*DSWAV 

where VREL = V + XAVG*R - ( (ZS-DSWAV) /2) *P and DELX is 
the length of each sidewall section, DSWAV is the draft of 
the ith sidewall section corrected for waves ,( ZS-DSWAV) /2 
is the average Z distance of the ith sidewall section from 
XREF where XREF is input as 9 feet forward of the transom. 
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These expressions are appropriate to the particular 
values associated with sidewall sections in either the port 
or starboard side. The yaw moment due to cross-flow drag 
is 

FND = FYD *XAVG 

where XAVG is the average X distance of the ith sidewall 
element form XREF. The roll moment is given by 

FKD = FYD * (ZS - DSWAV /2) 

The total force and moment values due to cross-flow drag 
are then 



FYD = 


FYD 


(1) 


+ 


FYD 


(2) 


FND = 


FND 


' (1) 


+ 


FND 


(2) 


FKD = 


FKD 


(1) 


+ 


FKD 


(2) 



where the subscripts (2) and (1) indicate the port and.- star 
board sidewalls respectively. 

The cross-flow drag coefficient is the value cor- 
responding to that for flow normal to a long flat plate, 
with a value of CDC = 1.28. 

2 . Slender Body Theory Terms 

Following the procedures indicated in reference 9 , 
together with evaluation of all integral terms involving 
derivatives of sectional properties by means of integration 
by parts the following results are found for a single side- 
wall (port or starboard, using appropriate values for each) 
The X-direction sidewall drag force is given by 
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FXH=-RH0/2*CDT * WAREA * U * U 



where WAREA is the wetted area of a sidewall section cor- 
rected for draft, roll, and pitch. CDT is the computed 
drag coefficient. 

The Z-direction sidewall force corrected for added 
mass effect is given by 

FZH = -G *BCO -U*U*A33S*TIIETA - U *A33S *W + Q * 

U * (A33S * XSS - BC2) - U*A33S * P *YLSW 

where BCO is the mass of the water displaced by the ith 
sidewall- element which causes a negative buoyant force, 

A33S is the vertical added mass of the ith sidewall element 
given by A33S = RHO * TT * BS **2/8 , where BS is its beam. 

BC2 is the total sidewall vertical added mass and YLSW is 
the lateral distance of the ith sidewall element form the 
centerline . 

The hydrodynamic lateral force is given by 

FYH = -A22S * U * (V + XSS * R - ZS * P ) 
where A22S is the lateral added mass at the stern given by 
A22S = RHO* .4 * PI * DSS**2/2 

where DSS is the sidewall draft. The moment about the Y-axis 
is given by 

FMH = -U *XSS*XSS* ASS3 * Q + G*BCOO + U * 

(A33 * XSS+BC2) * (W+U * THETA+ YLSW* P ) 
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where BCOO is associated with the sidewall volume between 
XREF and XS . This computation is available by entering a 
table of values computed offline in program SIDETAB and 
stored in common block WAVTAB by the INCON subroutine. 

The moment about the Z axis is given by 

FNH = FYH *XSS -U * ( (V-ZS *P) * BC5 + BC 6*R) 

where BC5 is the total lateral added mass and BC 6 is the 
moment associated with the total lateral added mass. 

The total slender body theory and crossflow forces, 
and moments for both sidewalls are then given by 



FX 


= 


FXH 


(1) 


+ 


FXH 


(2) 




FY 


= 


FYH 


(1) 


+ 


FYH ( 2 ) - FYD 




FZ 


= 


FZII 


(D- 


+ 


FZH 


(2) 




FK 


r= 


(FZH 


(2) 


- 


FZH 


(1) ) * YSW 


- ZS * FY + FKD 


FM 


= 


(FMH 


(1) 


+ 


FMH 


(2) ZS * 


FX 


FN 


=3 


(FXII 


(1) 


- 


FXH 


(2) ) * YSW + 


FNH (1) + FNH (2) +FND 



In the calculation of the individual sidewall force _ 
terms it is understood that the computation is appropriate 
to the instantaneous values of the separate elements as the 
vehicle moves . 
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I. SUBROUTINE STERN SEAL 

Since the bow and stern seals of the XR-3 are similar 
except for location the subroutines also are similar except 
for variable names . 

Reference to Figure 13 and the list below will acquaint 
the reader with the important variable names. 



Dimension 


Origin 


XS 


INCON ROUTINE BLOCK 225. XS is calculated 
by locating the CG and there calculating 
the distance to the transom. 


XSSI 


INCON ROUTINE BLOCK 500. This is input as 
the distance from the transom to the hinge 
point of the seal. 


ZS 


INCON ROUTINE BLOCK 200. ZS is calculated 
in the same manner as XS except that it 
measures the vertical distance from CG 
to the keel . 


ZSSI • 


INCON ROUTINE BLOCK 500. This is the ver- 
tical distance from the hinge to the keel. 


XSS 


This is calculated in INCON BLOCK 1302, 
using the information XSSI and XS and 
represents the horizontal distance from 
the center of gravity to the seal hinge. 


ZSS 


This is calculated in INCON BLOCK 1302 using 
the information provided by ZSSI and ZS 
and represents the distance from the center 
of gravity to the hinge. 


THSSI 


This is the angle of the seal face made 
with the craft vertical and is input in 
block 500 of INCON. 



Lift and drag forces in the stern seal result from the 
pressure differential between the seal and the plenum and 
the amount of the seal in contact with the water respectively. 
Since the plenum exhausts mainly under the stern seal we 
assume the bubble to extend to the tip of the stern seal. 
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Thus the stern seal becomes an air lubricated surface and 
the plenum pressure is extended all the way under the seal 
to the atmosphere. 

We can then assume that the forces acting upward con- 
centrate themselves at the point (B) in Figure 13 which re- 
presents the mid-point of the wetted surface. It is also 
logical to exert the drag forces at this same point, since 
they act along the wetted surface, the surface of friction. 

We find that the forces in the Z direction are much 
less than in the bow seal because the differential pressure 
is so much less than in the bow seal. The bow seal is acting 
against atmospheric pressure while the stern seal is acting 
against the plenum pressure. 

J. COMMON BLOCK TIE MAP 

The Common Block Tie Map shows which common blocks 
appear in each subroutine. 

The map is read by finding the common block in the left 
hand column and then reading to the right. An X appears in 
each column under the subroutine in which that common block 
appears . 

Some common ties have only one entry. This results 
from the removal of the SAM subroutine. 
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COMMON BLOCK TIE MAP 



OUtSKUU i JLiNi^ 

NAME 

COMMON 
BLOCK NAME 


MAIN 


AEROD 


BOWSL 


COLFIL 


FAN 


INCON 


INTGRL 


PROP 


RHS 


RUDDR 


SIDEWL 


STNSL 


WAVES 


AIR 


X 




X 




X 


X 






X 




X 


X 




AXIS 








X 




X 
















BMCO 


X 










X 


X 




X 




X 




X 


COLUMN 








X 




X 






X 










CONST 


X 




X 






X 




X 


X 


X 


X 


X 


X 


CNTRL 












X 






X 










CURVE 








X 




X 
















ENGINE 


X 










X 




X 


X 










EQNCO 


X 










X 


X 














FAERO 




X . 














X 










FAIR 




X 








XX 
















FANMAP 










X 


X 






X 










FORBS 






X 












X 










FORSS 


















X 






X 




FPROP 


X 














X 


X 










FROUDE 


X 










X 






X 










FRUD 


















X 


X 








GBOW 












X 






X 










GEOM 






X 






X 






X 




X 


X 


X 


GEOMBS 






X 












X 








X 


GF.OMSS 


















X 






X 


X 


GEOMSW 












X 










X 
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COMMON BLOCK TIE MAP (CONT.) 



bUBKUUTJ-Wii 

NAME 

COMMON 
BLOCK NAME 


MAIN 


AEROD 


BOWSL 


COLFIL 


FAN 


INCON 


INTGRL 


PROP 


RHS 


RUDDR 


SIDWL 


1 STNSL 

t 


WAVES 


GRAF 








X 




X 
















HEADG 








X 




X 
















KSWITCH 














X 




X 




X 


X 




LEAKER 






X 






X 












X 




MASSES 






X 






X 


X 




X 


X 


X 


X 


X 


MATRIX 












X 






X 










MSIDW 


















X 




X 






MWAVE 


















X 








X 


OPTION 












X 






X 










PLENUM 












X 






X 




X 




X 


PRIME 


X 










X 


X 




X 




X 






PRTINT 


X 


X 


X 




X 


X 


X 


X 


X 


X 


X 


X 


X 


PWAVE 












X 






X 










RISER ' 












X 














X 


ROLL 


X 










X 
















RUDDR 


X 










X 




X 


X 


X 








SIDE 












X 






X 




X 




X 


SOFTBS 




X 






X 


X 






X 










SOFTSS 










X 


X 






X 






X 




STABLE 












X 


X 














STSLR 












X 












X 




SUM 








X 




X 
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COMMON BLOCK TIE MAP (CONT.) 



SUBROUTINE 



NAME 

COMMON 
BLOCK NAME 


MAIN 


AEROD 


BOWSL 


COLFIL 


(fan 

1 


INCON 


INTGRL 


PROP 


RHS 


RUDDR 


SIDWL 


STNSL 


WAVES 


VALOLD 


X 










X 


X 




X 






X 




VARBLE 


X 


X 


X 




X 


X 


X 


X 


X 


X 


X 


X 


X 


WAVE 


X 




X 






X 






X 




X 


X 


X 


WAVEF 












X 
















WAVTAB 












X 










X 




X 



69 



IV . PROCEDURES 
’ — — — — ■ — 

A. OBTAINING STEADY STATE 

Following program modification to reflect the geometry 
of the XR-3, program runs were made to obtain steady state 
conditions. This was accomplished by comparing steady state 
program outputs of thrust and pitch angle with measured 
data obtained from the XR-3 testcraft. The testcraft mea- 
sured data consisted of steady state thrust and pitch angle 
measurements made at several velocities between 10 and 25 
knots in calm water. These measured data formed the basis 
upon which program modifications were to be made. 

It was first necessary to change the thrust versus 
velocity curve produced by the program so that it coincided 
with the measured data.- This was accomplished by adjusting 
the bubble frictional drag coefficient to bring the lower 
portion of the curves into agreement. A similar agreement 
was obtained at the high end of the curve by adjusting the 
skin frictional drag. This drag is a direct result of craft 
draft, and is therefore not easily obtained. Adjusting draft 
requires that both plenum pressure and fan flow be maintained 
in such a manner that the appropriate lift is attained 
without exceeding the horse power limitation of the fan 
engines. Successive runs using the constant speed, trim 
option of the program were made at various drafts until the 
appropriate values of thrust for a velocity was obtained. 



70 



Having found that draft which produces the required 
drag, hence, the proper thrust, bubble pressure and stern 
seal leakage were adjusted to bring the forces in the Z 
direction into balance. To accomplish this the plenum 
pressure that produced the required lift was found and then 
the stern seal leakage was adjusted to provide the proper 
fan flow which maintained fan engine horsepower within its 
limits. A small portion of the total seal leakage was 
assumed to exist at the joints between the bow seal and 
sidewalls, the remainder being under the stern seal. 

No attempt has been made to provide agreement between 
the craft transient behavior and that of the program. How- 
ever, the data herein presented appears to show that there 
is a fair amount of agreement in this area. The transient 
behavior of the craft and the simulation program provide 
a very fertile area for further investigation. 

B. DIFFERENCES ENCOUNTERED IN THE MODEL 

1 . Seal Configuration 

For reasons previously stated, the mathematical 
model of the 100B stern seal was selected to represent both 
seals on the XR-3. The forces ^calculated by the simulation 
program employing this model appear to be reasonable. That 
is, they are of the same order of magnitude as those scaled 
from the 100B. Leakage coefficients were not changed since 
the type of orifices employed on the XR-3 were considered 
to be similar to the 100B. 
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2. Plenum Dimensions 



For a first try at implementing craft plenum geo- 
metry actual craft measurements we re taken and used as input 
data to the simulation program. However, it was found that 
these dimensions did not produce appropriate XR-3 force 
and moment relationships. Through investigation of XR-3 
model data and subsequent program utilization the arbitrarily 
chosen nominal values of 10 and 20 feet, plenum width and 
length respectively, provided the desired results. Lift due 
to plenum pressure was highly sensitive to these plenum di- 
mensions. Therefore, when more accurate plenum pressure 
data becomes available these dimensions may require some 
refinement . 

C. SENSITIVITY TO CHANGES TO INPUT DATA 

It was found that the XR-3 is extremely sensitive to 
the initial conditions included in the input data set. Small 
disagreement between the required initial conditions and 
those input, produce violent transients that can cause the 
program to stop execution or to give erroneous results. 

Two of the most critical input parameters are craft 
draft and plenum pressure. This is true because the plenum 
pressure and draft are closely interrelated and cannot be 
separated for analysis. The plenum, having a large surface 
area, exerts the greatest single force on craft dynamics. 

Hence, a very small change in plenum pressure operating on 
this large surface area, 200 square feet, produces a significant 
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change in the Z direction force. Thus, imbalance of either 
draft or pressure, or draft and pressure, cause transients 
too great for the simulation to accept. These results also 
suggest that rapid maneuvers or inputs approaching step- 
functions will also produce transients too great to be ade- 
quately computed. 

It is recommended that an investigation be conducted 
to determine methods of improving the integration procedure 
to allow very rapid maneuvers. 

D. POSSIBLE ERRORS IN THE PROGRAM 

Many coefficients herein described were assumed to be 
correct. However, this may not be the case. Those coefv 
ficients that are calculated by their definitions such as 
Reynolds number and Froude number are most probable correct. 
However, leakage orifice coefficients, skin friction coef- 
ficients and the bubble drag coefficients have not been 
verified, though they seem to give adequate results. 

Other coefficients unknown to be exact are the moments 
of inertia. The simulation program provides a method for 

accurately calculating the moments of inertia, I I I 

X , Y f L 

and I . This option was not employed because accurate 

weight distribution data was not available. The moments of 

inertia I , I„ and I were obtained from references des- 
X Y z 

cribing the craft prior to many major modifications. Not 
the least of which was the installation of the auxiliary 
power supply and the replacement of bow and stern seals. 
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These modifications added a considerable percentage to the 

total craft weight. These weight additions were concentrated 

in several locations. Hence, the scaled values of 1^, , I y and 

I used in the simulation cannot be relied upon as being 
h 

accurate . 

No available data could be found concerning the product 

moment of inertia I . This moment of inertia is required 

X h 

when attempting to execute maneuvers involving turns or off 
centerline craft motions. Therefore, an estimated value 
using size scaling was employed and rendered reasonable craft 
motions under the above conditions. 

To determine the correct moments of inertia, it is re- 
commended that a thorough craft mass distribution analysis 
be conducted and the resulting parameters be input to the 
program utilizing the option mentioned above. 

A third area open to further improvement is that of 
refining the dimensions input to the sidewall map of the 
side geometry subroutine. The present representation for 
the XR-3 sidewalls loses much computational accuracy due to 
the large vertical stepsize presently used. 

A remapping of the sidewall geometry is required to 
more accurately describe the curvature in the lower regions 
of the sidewalls. This would result in the subroutine pro- 
viding more refined values of sidewall volume and surface 
area for the various XR-3 drafts and attitudes which are 
important for the calculation of sidewall -forces and moments. 
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Notwithstanding the above 
it is felt that the craft data 
that the simulation program is 
XR-3 motion. 



areas in which error exists 
and figures tend to indicate 
a valid representation of 
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V. PROGRAM VERIFICATION 



A. INSTRUMENTATION AND RECORDED DATA 

Instrumentation as presently installed in the XR-3 

V 

includes the following: 

1. Two strain gauges mounted directly on the engine 
mounts measure the thrust of the propulsion 
system. 

2. A rate Gyro mounted near C.G. is used to measure 
pitch and roll angles. 

3. A turbine type flowmeter which extends forward 
in undisturbed water measures velocity. 

Additional instrumentation presently operational on 
the XR-3 measures yaw, pressure in the plenum and seals, 
and rudder angles. The. output from these sensing devices 
are recorded on a 14 track Pemco Model 12 OB Magnetic tape 
recorder. Hard copies of the measurement data are obtained 
using a Hewlitt Packard two channel chart recorder. No 
filtering other than that natural filtering caused by the 
frequency response of the recorders was used with the re- 
sult that the data was quite noisy. The noise was consi- 
dered to be white and gaussian, therefore the mean value of 
the resultant curve was taken as the best estimate of the 
variable under cosideration . 

Shown in Table I are the data collected from test runs 
made by the XR-3 at Lake San Antonio, California, over a 
one-day period in calm water and no wind conditions. Each 
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entry represents a steady state condition. 



77 



Table I 



Test Craft Measured Data 



Velocity in kts. 


Thrust in 


8.5 


405 


9.0 


425 


10.5 


520 


12.0 


475 


13.0 


453 


14.0 


490 


14.5 


400 


15.0 


465 


16.5 


407 


17.0 


395 


19.0 


415 


19.5 


425 


20.0 


440 


21.0 


455 


22.0 


460 


22.5 


455 


23.0 


480 


24.0 


496 


24.5 


510 


25.0 


509 



pounds Pitch angle in degrees 

Not available 
Not available 
1.6 
1.55 

Not available 
1.4 

Not available 

. 8 

Not available 
Not available 
Not available 
Not available 

.2 

Not available 
Not available 
Not available 
Not available 
Not available 
0.1 

Not available 



78 



B. L&M PROGRAM DEVELOPED DATA 



1 . Steady State Trim Conditions 

Steady state trim conditions for the XR-3 using 
the L&M program were difficult to obtain. It was observed 
that the program initial conditions must be very close to 
the steady state values or the imbalance in forces generated 
within the program is so great that one of two cases occurs. 

Case 1 using Variable Time Step (R-K-M) 

The force imbalance causes such a radical transient 
that, in order not to lose curve information, the inte- 
gration step size is reduced until the minimum allowed step 
size is reached and the program execution stops. 

Case 2 using Fixed Time Step (R-K) 

The danger in this method is that the step size may be 
too great (note: in some cases, .01 seconds is too large) 

and curve information is lost. When information is lost 
the program calculations begin generating erratic number 
values thus causing an underflow or overflow condition to 
exist within the program which stops program execution. Thus, 
care must be exercised in selecting a stepsize to be used. 

The results shown in Table II were obtained using the fixed 
step size option with a step size of .001. The constant speed 
and variable thrust options of Block 1 of the INCON Subroutine 
were selected. Since the force equations of the program are 
dependent upon speed, using these options removes the surge 
equation from the computations and adjusts thrust to balance 
the forces in the X direction. 
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Table II. L&M Program Initial Conditions 



Speed 
(in knots) 


Thrust 
(in pounds) 


Pitch 7\ngle 
(in degrees) 


Draft 

(in inches) 


10.0 


504.6 


1.72 


6.6 


12.5 


419.3 


1.31 


6.5 


15.0 


392.4 


. 94 


6.35 


17.5 


397.8 


. 6 


6.15 


20.0 


424.7 


. 29 


5.98 


22.5 


464.7 


.13 


5.78 


25.0 


512.8 


.09 


5.56 


27.5 


567.3 


.06 


5.3 


30.0 


625.5 


. 05 


3.02 



Plenum pressure 24.86 psf 
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To balance the forces in the Z direction short time 
runs of .01 seconds or less were used. For each run ad- 
justments were made to the initial draft and/or plenum 
pressure to bring these forces into balance. Care was 
taken to avoid fan flows and pressures above that level 
requiring greater than 15 horsepower, which is the maximum 
horsepower available from the fan engines. Continuing in 
this manner the initial trim conditions of draft, pitch 
angle, and thrust required were obtained for each of the 
speeds shown in Table II. Once these initial conditions 
were known, the step size was increased to .01 and it was 
found that no curve information was being lost indicating 
that the transient was of a less violent nature. Long time 
runs using the variable time step size verified that these 
values were' the correct initial trim conditions. A straight 
line segment was used to join the points so that linear in- 
terpolations may be used between the tabular values. 

To obtain good agreement between measured and simulated 
data required the manipulation of two quantities that pro- 
duce total ship drag; the bubble drag, FXPWAV, and the skin 
friction drag FX. Adjustment of the forces in the X direction 
was accomplished by adjusting draft to provide the proper 
drag force at the upper end of the velocity curve where the 
frictional drag component exerts the greatest influence over 
total ship drag. The next step was to adjust the coefficient 
of the bubble drag to provide the minimum point in the total 
drag curve at the velocity indicated by Table II , about 17 . 0 
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knots, and to also provide the required drag at the low end 
of the drag curve where the bubble drag exerts the greatest 
influence . 

To obtain agreement between measured and simulated 
pitch angle required some programming changes. Reference 11 
indicates that a pressure gradient exists in the plenum 
chamber. This gradient is such that as the craft increases 
in speed the center of pressure moves aft. This causes the 
craft to have a pitch up attitude at slow speeds and a gra- 
dually smaller positive pitch as speed increases. Using the 
experimentally measured data shown in Table I , a set of points 
representing the craft center of pressure position as a 
function of speed was calculated which would produce the 
appropriate craft pitch angle for a given speed. A linear 
function used to approximate these points was then included 
in the simulation program. 

Figures 14, 15 and 16 show plotted points of the two 
tables of 'data. Triangles represent experimentally deter- 
mined data and squares the program output. Clearly the pro- 
gram simulates quite well the thrust and pitch angle curves. 
Accurate data on draft was not available from XR-3 tests 
due to a defective height sensing device. However, it 
is felt that the draft data when available will closely 
fit the curve shown. 

2 . Presentation of Unverified Program Data 

Figures 17, 18 and 19 demostrates the use of the 
variable thrust option. Thrust is input to Subroutine PROP 
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Figure 14 
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Figure 16 
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Figure 17 
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PITCH ANGLE (x5.0E 02 degrees) 



Figure 18 

THRUST CHANGE: Pitch Angle Versus Time 
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DISPLACEMENT (inches) 
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Figure 19 

THRUST CHANGE: Z Displacement versus Time 
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utilizing the Block 16 variable thrust option of Subroutine 
INCON. This thrust is applied equally to both screws and 
is mapped as follows. Thrust for 20 knots is initially 
applied and maintained for 10 seconds. Between time 10 
seconds and time 15 seconds the craft experiences a ramp 
function of thrust building up to the tabular value for 
25 knots. This thrust is maintained for the remainder of 
the run. 

Figure 18 shows that the pitch angle follows the ex- 
pected curve, in that initially pitch up is experienced 
as additional thrust is applied. This is followed by gra- 
dually dampening oscillations toward the value of pitch 
angle expected at 25 knots. 

Figures 17 and 19 show that speed through the water 
increases and that draft decreases, as expected. Addi- 
tionally, examination of Block 16 of Subroutine INCON makes 
it clear -that the engines may be unevenly thrusted rather 
than equally thrusted as shown. 

Figures 20, 21, 22 and 23 show the behavior of the 
program for a turning maneuver. The rudder input vs. time 
is shown in Figure 20. Corresponding roll angle vs. time, 
pitch angle vs. time, and the ship horizontal plane motion 
are in Figures 21, 22 and 23 respectively. 

It is interesting to note that the roll and pitch 
damping coefficients appear to be different during the time 
the rudder is off centerline than when the rudder is re- 
turned to centerline. However, these differences are to 
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20 Knot Turn: Rudder Angle vs. Time 
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Roll Angle (x5.0E-01 degrees) 

-020 -025 ' -020 -015 -010 -005 000 COS 



Figure 21 




Fig. 20 Knot Turn: Roll Angle vs. Time 
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Pitch Angle (xl.OE-Ql degrees) 
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Figure 22 
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Displacement (x2.0E+02 feet) 



Figure 23 




X Displacement (x2.0E+02 feet) 

20 Knot Turn: X Y Displacement 
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be expected due to the nonlinearity of the system. The 
oscillation appears to be about 0.5 H^. This oscillation 
was not felt by the authors during craft test runs. However, 
it may not be sensed due to its high frequency and small 
amplitude. Unfortunately at this writing no accurate turn 
data is available. 
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VI. DISCUSSION AND EVALUATION OF SIMULATION RESULTS 



A. DISCUSSION OF PRODUCT MOMENT OF INERTIA; I v _ 

X h 

No information has been ^available on the value of the 

product moment of inertia I v „ for the XR-3. Therefore, 

throughout the XR-3 program development the moment of inertia 

I was set equal to zero. This decision had no effect on 
X U 

craft motions for calm water straight runs since craft 
motion was restricted to translational movement in the Z 
and X direction and rotation about the Y axis. However, 
when turns were attempted and waves off the bow were intro- 
duced, program instabilities were observed. These may be 
seen by referring to Figures 24, 25 and 26 which show that 
for small perterbations the craft reacts in a stable manner, 
but as thd rudder angle increases toward 20 degrees the craft 
becomes unstable even though waves have not yet been intro- 
duced. A weight scaled value of I was obtained from the 
100B input data set and was later used as a representative 
value for the XR-3. 

B . 5 ' DEGREE TURNS WITH AND WITHOUT WAVES 

In an attempt to verify conditions reported in Reference 
8 for the 100B model, a 20 knot run with a 5 degree rudder 
turn with and without waves was introduced. The wavelength 
and wave height were scaled from Reference 8 according to 
relative craft size. Hence, a wavelength of 72 feet and a 
wave height of 0.16 feet were used with the XR-3 model. Wave 
amplitude and frequency of encounter are shown in Figure 29. 
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■ ROLL ANGLE (degrees) 



Figure 24 

FIVE , FIFTEEN AND TWENTY DEGREE TURNS: 
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PITCH ANGLE (degrees) 



Figure 25 

FIVE, FIFTEEN AND 20 KNOT TURNS: 
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YAW ANGLE (x5.0E 01 degrees) 



Figure 26 



FIVE, FIFTEEN AND TWENTY DEGREE TURNS: 
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It is interesting to note that this same frequency can be 
seen in roll angle and pitch angle as shown in Figures 27 
and 28. It is observed in Figure 27 that the amplitude of 
roll increases almost linearly as the wave front moves away 
from the ship's head, as would be expected for beam seas. 

C. TURNS WITH 5, 15 AND 20 DEGREE RUDDER 

Figures 31, 32, and 33 show the roll, pitch, and yaw 

motions produced by 5, 15, and 20 degree rudder turns at 

20 knots in calm water with the product moment of inertia 

1^2 inserted. In all cases, damping and system stability 

were observed. Similar turns made with I equal to zero 

X z 

revealed the presence of roll and pitch instability for 
large rudder angles. In both cases where 5 degx*ee rudder 
turns were made the variations in the i~oll and pitch angles 
remained small and decayed to zero within 15 seconds or 
about 5 cycles. 

D. COMPARISON OF SIMULATION AND MEASURED DATA 

At the time of this writing measured pitch and roll 
data was not in a form suitable for definitive comparisons. 
However, recorded data of testcraft motion during turn ma- 
neuvers as shown in Figure 34 tend to validate the presence 
of damped roll and pitch oscillations which are also pre- 
sent in the simulated turns as shown in Figures 31, 32 and 33. 

Measured and recorded plenum pressures tend to cor- 
respond with those calculated by the simulation program. 

This pressure appears to vary between 24-25 pounds per square 
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ROLL ANGLE (xl.OE 00 degrees) 



Figure 27 

FIVE DEGREE TURN: Roll Angle Versus Time 
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PITCH ANGLE (x2.0E-01 degrees) 



Figure 28 

FIVE DEGREE TURN: Pitch Angle versus 

Time with and without waves 
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•WAVE HEIGHT (x5.0 E-02 feet) 



Figure 29 

FIVE DEGREE TURNS: Wave Height versus Time 



1. Turn with waves 

2. Turn without waves 
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S PLACEMENT (x5.0E 01 feet) 



Figure 30 



FIVE DEGREE TURN: X Displacement Versus Y Displacement 
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ROLL ANGLE (deg 



Figure 31 

FIVE, FIFTEEN AND TWENTY DEGREE TURNS: 
Roll Angle versus Time 



1. 20 degree rudder 

2. 15 degree rudder 

3. 5 degree rudder 



n 

c-> 



o 

o 

o 



] v 



if) 

Q) 

O 

M 



• -t 

o 

o 

i 



tv 

O 

cj 

i 



o* 

o 

o 



to 

o 

o 

» 



TIME (seconds) 




104 



PITCH ANGLE (xl.OE-Ol degrees) 



Figure 32 

FIVE, FIFTEEN AND TWENTY DEGREE TURNS: 
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YAW ANGLE (x5.0E01 degrees) 



Figure 33 

FIVE, FIFTEEN AND TWENTY DEGREE TURNS : 
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foot throughout all the various monitored craft motions. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



The conversion of the L&M program from the 10 OB model 
to the smaller XR-3 model involved considerably more time 
and effort than was anticipated at the onset. A great deal 
of time was required for familiarization with the program, 
for providing a more useful output format and for adding 
plotting routines. Considerable effort was required to 
obtain the straight run, calm water trim conditions for 
various speeds. Plenum center of pressure movement with 
velocity was one of the more difficult problems encountered 
in obtaining the proper initial trim conditions. It is 
recommended that future investigators write an initial 
condition program to assist in reducing the time and effort 
required to arrive at the proper trim conditions. 

Good agreement between measured and computed data for 
calm water straight runs was possible by the adjustment of 
the bubble drag coefficients, draft, bubble pressure and 
stern seal leakage. The measurement data for turning ma- 
neuvers could not be used for a check on the simulation in 
the form in which it is presently recorded. It is recom- 
mended that the recorded data be processed digitally and 
suitably reduced to a form useful for comparison with the 
L&M program output. Qualitative comparison of the XR-3 turn 
maneuver data however, has indicated that pitch and roll 
motion have a lightly damped response quite similar to the 
L&M program output. Much work remains to be done in checking 
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the accuracy of the computed transient response. 

Additional qualitative checks were made on the com- 
puted XR-3 turn maneuvers by comparison with the 100B 
model response. The results indicated similar transient 
characteristics in which damping decreases as pitch or 
roll angle increase. 

The simulation study of the XR-3 has demostrated the 

importance of the product moment of inertia I in main- 

X z 

taining stability in roll and pitch for rudder angles 

greater than 5 degrees. Future investigators should attempt 

to obtain an accurate value for I . This can be done by 

x z d 

utilizing the discrete mass distribution option in the 
INCON Subroutine. 
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APPENDIX A 

SES MOTIONS AND LOADS PROGRAM 
USERS MANUAL 
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SES MOTIONS AND LOADS PROGRAM 
USERS MANUAL 
by 

James Bentson, Theodore P. Sargent and Alfred 
Oceanics Incorporated 
with modifications by 

LCDR Don G. Leo and LT Richard Boncal U. S. 



I. Raff 



Navy 
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INTRODUCTION 



The SES Motions and Loads Program is a Fortran IV pro- 
gram capable of calculating the time histories of the motions 
and loads of a surface effect ship operating in the on bubble 
mode in a seaway in six degrees of freedom. It is programmed 
in a "modular" fashion, i.e. each major element of the craft 
has its own subprogram. The input is in the form of basic 
data, i.e. sidewall offsets, seal geometry, fan map data 
points, etc. and enables the user to make a wide range of 
parametric changes to the particular craft being modeled 
by changing only the input data. The program has the ca- 
pability of printing and/or plotting a wide range of var- 
iables as a function of time; such as forces and moments 
from each element, total craft forces, moments and ac- 
celerations, leakage gaps and flow rates, fan flow rates 
and power required, etc. The current program has been 
checked out and run on an IBM 360 computer (using version 
Release 20.6.) at .the W. R. Church Computer Facility, NPS , 
Monterey, Ca. 

This report is intended to be primarily a users manual 
for the program as it exists. The details of how one would 
rewrite various modules to reflect changes in the craft de- 
sign are discussed extensively in Appendix A-A. 
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PROGRAM ORGANIZATION 



Within a given run the program proceeds as follows. 
First the main program calls the input and initialization 
subroutine INCON. The program then proceeds to calculate 
the motion time histories by calling the various modular 
subroutines and the integration routine for as many time 
steps as is necessary to complete the run. During the 
course of the calculation the various subroutines will 
output those variables selected by the print option switches 
at the appropriate input print interval. In addition, 
data may be written on scratch files to be used for sum- 
marizing the output and doing bending moments and shear 
calculations. At the present time the Shears and Moments 
Subroutine has been removed, however, it may be reintroduced 
at any time with a minimum of changes. 

When the run is completed, the program, depending on 
the options chosen, may then print and/or plot the output 
summary and proceed to do the bending moment and shear cal- 
culations. After these steps are completed the main pro- 
gram returns to INCON to read data for the next case. 

The above procedure is repeated until all cases are 
finished. . 
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PROGRAM REQUIREMENTS 



As run on an IBM 360 computer under OS/360MVT version 
(Release 20.6) , the program requirements are as follows; 

TOTAL CORE TIME 

Program length 201,456 -^q NA 

Core size required for execution 250,000-. 0 SEE 

BELOW 



The total core required for a given step includes the 
program, system routines, data storage, etc. The time re- 
quired to run a single case depends on the nature of the 
case and the length of time to be simulated. The ratio of 
required execution time to real time has been found to be 
as low as 1/3 for some calm water cases. Higher ratios 
occur where the transient curves are rather sharp. 

The program currently uses the following FORTRAN Data 
Set Reference Numbers (Unit Numbers) : 



Unit Number 
1 
2 

5 

6 

10-15 



Usage 

Scratch file for COLFIL plotting package. 

Scratch file for lateral plane output 

summary 

Input 

Output 

Sidewall Integral Table data (see Appendix A--C) 
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IBMF 



Intermediate file containing time histories 
required by bending moment calculations. The 
unit number IBMF is read in using the 01501 
input card (see discussion of Input) 

INPUT AND OUTPUT 

The input for the program is handled by subroutine INCON 
and is discussed there. The output is controlled by the 
individual subprograms and is discussed separately there. 
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DESCRIPTION OF SUBPROGRAMS 
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Main Program 



Discussion - The main program contains the logic for inter- 
connecting the various subroutines. In addition, it com- 
pares the running value of time with the finish time, 
calculates the next value of time for printing after each 
print time, does some output, etc. The main program also 
contains the logic for the integration of the trajectory 
(x, y,MM by the trapazoidal rule. 

Input - None 

Output - The output from the main program is controlled by 
the print option switch ITRAJ. If the value of this switch 
is 1, the program will print the current values of the time 
and the craft velocities and displacements (angular and 
translational) as shown in the sample output at the end of 
the report . 
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Subroutine AEROD 



Discussion - Subroutine AEROD calculates the aerodynamic 
forces and moments on the craft. The subroutine currently 
uses curve fits of the aerodynamic coefficients vs sideslip 
angle. Appendix AhA discusses how the program may be changed 
to input tabular data for the aerodynamic coefficients. 

Usage - CALL AEROD 

Input - The input data for AEROD is controlled by block 10 
of subroutine INCON and is discussed there. 

Output - The output from AEROD is controlled by print option 
switch IAEROD . If this switch is 1, the program will print 
the total forces and moments calculated as shown in the 
sample output at the end of the report. 
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Subroutine BOWSL 



Discussion - Subroutine BOWSL calcualtes the forces and 
moments due to the bow seal and the leakage flow rate 
associated with any gaps which open up under the various 
seal stations. 

Usage - CALL BOWSL 

Input - The input data for BOWSL is controlled by blocks 
3, 6, and 17 of subroutine INCON and is discussed there. 

Output - The output from BOWSL is controlled by print option 
switch IBOWSL. If this switch is 1, the program will print 
the gaps and wetted lengths at the various stations on the 
bowseal as well as the total forces and moments due to the 
bowseal. A sample print point is given in the sample out- 
put at the end of the report. 
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Subroutine COLFIL 



Discussion - Subroutine COLFIL is used to format the pro- 
gram output. It can provide up to 16 variables listed in 
two tabular summaries of 8 variables each. The subroutine 
has a plotting package that permits use of either the print 
plot or calcomp (continuous line) option. Multiple curves 
may be plotted and up to 10 graphs may be drawn on each run. 
Print switches control the printing by each subroutine of 
its forces and moments. This option is very useful in de- 
bugging the program. Normally it provides too much output 
and is too cumbersome. 

Usage - CALL COLFIL 

Input - The input is controlled by the control card with 
the control tag 00104 and is discussed in the description 
of Subroutine INCON in the block 1 control card section. 
Additional input required for the plotting package is con- 
tained in the Users Manual for the COLFIL routine. 

Output - For output format consult the COLFIL plotting 
package users manual . 
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Subroutine DMINV 



Discussion - Subroutine DMINV is a matrix inversion package 
used for inverting the mass matrix. The method used is the 
standard Gauss-Jordan technique. The determinant of the 
matrix is also calculated. If the determinant is zero the 
matrix is singular and the program vzill stop (see Appendix 
A-B Abnormal Termination Messages) . 

Usage - CALL DMINV (A, N, D) 

Description of Parameters 

A = The input matrix, destroyed in computation and 
replaced by resultant inverse 
N = The order of matrix (maximum = 6) 

D = The resultant determinant of the matrix 

Input - None 

Output - None 



124 



Subroutine FAN 



Discussion - Subroutine FAN is used to calculate the inflow 
to the plenum, stern seal, and bowseal from the fans for 
a given pressure differential across the fans. The sub- 
routine uses tabular data and allows for changes in fan 
speed as well as the number of fans. Fan speed changes are 
currently handled by using the fan scaling laws rather than 
reading data as a function of fan speed. 

Usage - CALL FAN 

Input - The input for subroutine FAN is controlled by block 
19 of subroutine INCON and is discussed there. 

Output - The output from FAN is controlled by print option 
switch IVEL. If this switch is 1, the program will print 
the flow rates and pressure differentials across the bow 
seal, cushion and stern seal fans, as well as the total fan 
horsepower required as shown in the sample output at the 
end of the report . 
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Function FG1 



Discussion - The purpose of FG1 is to perform the task of 
evaluating a tabular function y = f(x). The procedure used 
is linear interpolation between data points. If the input 
value of the independent variable does not fall in the range 
of its input tabular values the appropriate table extreme 
is held. The function y must be single valued and the ta- 
bular values for the independent variable must be stored in 
ascending order i.e. XTAB (J) , XTAB (J+l) . The increment be- 
tween tabular values need not be constant. 

Usage - Y = FG1 (X, N, XTAB , YTAB , IMEM) 

Description of Parameters 

Y = The output value of the tabular function 
X = The value of the independent variable to be used 
to evaluate the function y = f(x) 

N = The number of elements in the XTAB and YTAB arrays 
XTAB = The name of the array that contains the values of 
the independent variable in ascending order 
YTAB = The name of the array that contains the values of 
the dependent variable corresponding to the values 
of the independent variable stored in XTAB 
IMEM = Dummy variable used to store location of last entry 
into table, to be used as initial point of search at 
next entry 

Input - None 
/ 

Output - None 



126 



Subroutine FORIT 



Discussion - Subroutine FORIT is used to obtain a Fourier 
analysis of a periodically tabulated function. It computes 

the coefficients of the desired number of terms of the Fourier 

<-M 

series f (x) = a(0) + £- (a(k)cos kx + b(k)sin kx) , to approx- 

k=l 

imate the given set of periodically tabulated values of a 
function. 



Usage - CALL FORIT (FNT, N, M, A, B, IER) 

Description of Parameters 

FNT = Vector of tabulated function values of length 2N+1 
N = The interval such that 2N+1 points are taken over 

the interval (0,2TI). The spacing is thus 2TT/(2N+1). 
N must be greater than or equal to M 
M = Maximum order of harmonics to be fitted. M must 
be greater than or equal to zero 
A = Resultant vector of Fourier cosine coefficients 
B = Resultant vector of Fourier sine coefficients 
IER = Resultant error code where: 

IER=0 No error 
IER=1 N less than M 

IER=2 M less than 0 

Input - None 
Output - None 
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Subroutine INCON 



Discussion - Subroutine INCON contains the logic for the 
reading of all input data, the initialization of variables 
and for the initiation of new cases. The subroutine is pro- 
grammed so that the input data is non-structured i.e. there 
is no set order of input cards from one read statement to 
another. However, a given read statement may have more than 
one input card. In runs with multiple cases, only those 
input parameters which are to be changed for the new case 
need be read in. 

Usage - CALL INCON 

Input - The input logic is subdivided into twenty- two ’'block 
Each block being for a particular subroutine or function e.g 
stern seal, sidewall, etc. Switching to the proper input 
block is done by means of a control tag on the input cards. 
Once control has been transferred to a particular block ad- 
ditional input may be read by specific format statements if 
necessary, however, the card which contains the control tag 



may also be used 


to input up to 


seven variables for 


the 


selected block. 


The format of 


the input card used 


for swi- 


tching is as follows: 






Column 


Format 


Entry 




1-3 


13 


control tag 




4-5 


12 


option tag 




6-75 


7F10 


input variables 


(if any) 



The control tag used to transfer control to the appro- 
priate block is of the form OMN where M and N are integer 
constants, i.e., if the control tag is 015 the program will 
go to statement 1500. 

The option tag is loaded into a variable called IOPT 
and allows different branches of logic within a given input 
block. If no branches are used within the input block the 
option tag may be zero. 

The input variables are stored in an array called TEMP 
and may be set to a given parameter by using an arithmetic 
replacement statement in the appropriate block (see various 
blocks in attached listing for examples) . 

Example: When read, a card of the form 



Column 


3 


5 




15 


25 


35 


80 




ir 




• 




•l? 








017 


02 


3.5 


blank 




1.0 


blank 



will transfer control to statement 1700. Within statement 
1700 the logic might look as follows 
1700 CONTINUE 

GO TO (1710, 1720), IOPT 
1710 A = TEMP (1) 

B = TEMP (2) 

GO TO 10 

1720 C = TEMP (1) 

D = TEMP (2) 

E = TEMP (3) 

GO TO 10 
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In the above example the option switch IOPT is 02 and 
will cause the program to go to statement 1720 which will 



then set C=3.5, D=0 . 0 , E=1.0 and then will return to statement 
10 which is the read statement for the input cards. 

The various input blocks as currently programmed are as 



follows : 




Block Number 


Function 


1 


Executive input (time steps, print switches, 
option switches, etc.) 


2 


Mass and inertial input 


3 


Craft geometric input 


4 


Sidewall input 


5 


Stern seal input 


6 


Bow seal input 


7 


Plenum input 


8 


Propulsion geometry input 


9 • 


Rudder input 


10 


Aerodynamic input 


11 


Wave input 


12 


Initial conditions 


13 


End of case indicator 


14 


End of run indicator 


15 


Not used 


16 


Thrust, rudder motion input 


17 


Bow and stern seal pressure input 


18 


Title card 


19 


Fan map input 
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Card 

1 



Option 

Card 

1 

2 



Columns 

1-5 

6-15 

16-25 

26-35 

36-45 

46-55 

tag = 2 
Columns 
1-5 
5 

10 . 

15 

20 

25 

30 

35 

40 



Format 

I 

F 

F 

F 

F 

F 



Entry 

00101 - control tag and option 
Start time (generally zero), sec. 
Finish time, sec. 

Initial integration time step, sec. 
Print interval, sec. 

Print start time, sec. 



Print switches (if=l print if=0 no print) 

Format Entry 

I 00102 - control tag and option 

I IACCEL, print switch for lateral 

acceleration and integration time 
step 

I ■ IVEL, print switch for subroutine FAN 

I ITRAJ, print switch for main program 

I ISIDWL, print switch for subroutine 

SIDEWL 

I IBOWSL, print switch for subroutine 

BOWSL 

I ISTNSL, print switch for subroutine 

STNSL 

I I WAVES , print switch for subroutine 

WAVES 

I IRUD, print switch for subroutine 

RUDDER 
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Option 


tag = 2 


(Continued) 




Card 


Columns 


Format 


Entry 


2 

(Cont . ) 


45 


I 


IPROP , print switch for sub- 
routine PROP 




50 


I 


IAEROD, print switch for sub- 
routine AEROD 




55 


I 


IRHS , print switch for subroutine 
RHS 


Option 


tag = 3 


Integrator 


Tolerances 


Card 


Columns 


Format 


Entry 


1 


1-5 


I 


00103, control tag and option 


2 


1-2 


12 


Number of integrators to be used 
(maximum of 15) 




3-4 


12 


Control tag for integration pro- 
cedure 


3, 4. . . 


1-80 


8F10 


Integrator error tolerance, one 
per integrator 



Note: The usual procedure for selecting integrator tolerances 

_ r __ O 

is to choose values on the order of 10 3 to 10 . However, 

for each new craft the following procedure is used; make a 
run with a given set of tolerances; halve the tolerances and 
rerun. If the results do not change drastically, try doubling 
the original tolerances and running. As changes are noted in 
the integrator outputs (i.e. craft velocities, etc.), hold the 
tolerances for those integrators fixed and vary the others 
until the largest tolerances possible are obtained. Once 
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chosen for a given craft the tolerances rarely have to be 
changed. A control tag (JQQ in the program) is used to se- 
lect integration procedure. A blank provides the original 
R-K., procedure. A"one” provides a fixed step size runga- 
kutta procedure that could be useful in selecting integrator 
tolerances. A "two” provides runga-kutta variable step size, 
but does not print out the step varying procedure. Only the 
minimum stepsize is output. 



Option tag = 4 
Ca^d Colum ns 

1 1-5 

2 5 
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Tabular Summary of Output 
Format Entry 

I 00104, control tag and option 

I Print switch for summary one 

1 = print, 0 = no print 
• Print switch for summary two 

1 = print, 0 = no print 



Option tag = 5 
Card Columns 

1 1-5 

6-15 

16-25 



Program Option Switches 
Format Entry 

1 00105, control tag and option 

F Switch for lateral plane motions 

1 = only lateral plane motions 

0 = six degrees of freedom 

F Switch for constant speed 

1 = surge equation not used 

0 = surge velocity allowed to vary 
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Card 


Columns 


Format 


Entry 




26-35 


F 


Switch for trim 



1 = thrust varied to give proper 
speed. If this option is se- 
lected, then there must be input 
in Block 16 an estimate of thrust 
and a 0.0 for the number of data 
points . 

0 = thrust as input 



Block 2 


- Mass 


and Inertial 


Input 


Option 


tag = 1 


Summary Mass Properties Card 


Card 


Columns 


Format 


Entry 


1 


1-5 


i 


00201 - control tag 




6-15 


F ' 


Total craft weight, pounds 




16-25 


F 


Longitudinal center of gravity 
feet forward of transom 




26-35 


F 


Vertical center of gravity 
feet above baseline (keel) 




36-45 


F 


I v „ , mass moment of inertia 

2 

about x axis, slug - ft. 




46-55 


F 


I , mass moment of inertia 

2 

about y axis, slug - ft. 




56-65 


F 


I , mass moment of inertia 
zz 

2 

about z axis, slug - ft. 




66-75 


F . 


I , mass moment of inertia 
xz 

2 

about xz axis, slug - ft. 
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Note: If this option is used, shears and moments will not 

be calculated correctly. 

Option tag = 2 Discrete Mass Distribution Cards 



Card 


Columns 


Format 


Entry 


1 


1-5 


I 


00202 - control tag 


2, 3, 


4 . . . 1 -10 


F 


Discrete weight, pounds 




11-20 


F 


Longitudinal c.g., feet forward 




21-30 


F 


of transom 

Transverse c.g., feet to starboard 




31-40 


F 


Vertical c.g., feet above baseline 


Last 


1-10 


F 


(keel) 

-1.0 - control tag 


card 









The mass distribution is specified by discrete weights at 
various locations with one card used for each weight. Sym- 
metry about the longitudinal center-plane is assumed, so that 
only weights on the starboard side, with atransverse location 
greater than zero should be specified. Thus, the total weight 
on the input cards is one half the total craft weight. The 
last card is used to signal the end of the discrete weights. 

The maximum number of discrete weights is 201. 

Block 3 Craft Geometric Input 
No option tags 

Card Columns Format Entry 
1 1-5 I 00300 - control tag 

6-15 F Number of stations on port sidewall 
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Card 


Columns 


Format 


Entry 






16-25 


F 


Number of stations 


on starboard 








sidewall 






26-35 


F 


Number of stations 


on bow seal 




36-45 


F 


Number of stations 


on stern seal 




46-55 


F 


Total craft length 


, f t . 


Note : 


The maximum 


number 


of stations is 11 for 


each element. 



The program will automatically subdivide the element into 
N-l sections (N is the number of stations) and calculate the 
x and y coordinates of each station. 



Block 4 Sidewall Input 
Option tag = 1 Sidewall Input 



Card 


Columns 


Format 


1 


1-5. 


I 




6-15 


F 


1 


16-25 


F 




26-35 


F 




36-45 


F 




46-55 


F 



Entry 

00401 - control tag and option 
y distance from centerline to 
sidewall, ft. 

Average wetted length of side- 
wall, ft. 

Leakage orifice coefficient of 
sidewall 

Cross-flow drag coefficient of 
sidewall 

Average beam of sidewall, ft. 
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Note: Under this option tag the program will also read the 

sidewall integral table file stored on Unit 10. 



Block 5 Stern Seal Input 
No option tags 
Card Columns Format 
1 1-5 I 

1 6-15 F 



16-25 



26-35 F 

36—45 F 

46-55 F 

56-65 F 

66-75 F 



Entry 

00500 - control tag 
x coordinate of seal hinge, ft. 
forward of transom 
z coordinate of seal hinge, ft. 
above keel 

Base leakage area, sq. ft. 

Seal leakage orifice coefficient 
Angle between leading edge of seal 
and craft vertical, deg. 

Pressure differential between 

seal bag and bubble, psf 

Length of leading edge of seal, ft. 



137 



Block 6 Bow Seal Input 



No option tags 

Card Columns Format 



1 1-15 I 

6-15 F 



16-25 F 

26-35 F 



36-45 F 



46-55 F 



.56-65 F 



66-75 F 



Entry 

00600 - control tag 

x coordinate of seal hinge, ft. 
forward of transom 

Seal leakage orifice coefficient 

Pressure differential between seal 

bag and bubble, psf 

Z coordinate of the seal hinge 

above keel , feet 

Angle between leading edge of 

seal and craft vertical, deg. 
Length of the leading edge of the 
seal, feet 

Base leakage area, ft. 



Block 7 Plenum Input 

Option tag = 1 Plenum Geometry 
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Card 


Columns 


Format 


Entry 


1 


1-5 


I 


00701 - control tag and option 


1 


6-15 


F 


Plenum length at water surface, ft 




16-25 


F 


Plenum width at water surface, ft. 




26-35 


r- 


Not used at present 




36-45 


F 


Plenum width at deck, ft. 




46-55 


F 


Plenum length at deck, ft. 




56-65 


F 


x coordinate of center of pressure 
ft., forward of transom 




66-75 


F 


Plenum average height 


Option 


tag = 2 


Critical Froude Number 


Card 


Columns 


Format 


Entry 


1 


1-5 


I 


00702 - control tag and option 




6-15 


F 


Froude number corresponding to 
hump speed 


Block 


8 • Propul 


si on Input 




Option 


tag = 1 


Basic Data 




Card 


Columns 


Format 


Entry 


1 


1-15 


F 


00801 - control tag and option 




6-15 


F 


x coordinate of propeller center, 
ft. forward of transom 




16-25 


F 


y distance from centerline to 
propeller center, ft. 




26-35 


F 


z coordinate of propeller, ft. 



above keel 
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Block 


9 Rudder 


Input 




Option 


tag = 1 


Rudder Geometry 


Card 


Columns 


Format 


Entry 


1 


1-5 


I 


00901 - control tag and option 




6-15 


F 


x coordinate of centroid of rudder, 




16-25 


F 


ft. forward of transom 
y distance from centerline to rudder 




26-35 


F 


centroid, ft. 

z coordinate of centroid of rudder, 




36-45 


F 


ft . above keel 
Rudder span, ft. 




46-55 


F 


Rudder aspect ratio 




56-65 


F 


Rudder area, sq. ft. 




66-75 


F 


Average Thickness ratio of rudder 


Block 


10 Aerodynamic Input 


section 


No option, tags 
Card Columns 


Format 


Entry 


1 


1-5 


I 


01000 - control tag 




6-15 


F 


Reference length, ft. 




16-25 


F 


Reference width, ft. 



Note: The reference length and width are those used to 

nondemensionalize the wind tunnel data used in the program. 
Block 11 Waves Input 

Option tag = 1 Amplitude vs. frequency 
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Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01101 - control tag and option 




6-15 


F 


Number of wave components , max- 
imum of 10 




16-25 


F 


Initial wave heading relative to 
craft headings (180° = head seas' 
deg . 


2, 3, 4. 


. . 1-10 


F 


Wave frequency, rad. /sec. 




11-20 


F 


Wave amplitude, ft. 


Option 


tag = 2 


Amplitude 


vs. wave length 


Card 


Columns 


Format 


Entry 


1 


1-5 


I 


00102 - control tag and option 




6-15 


F 


Number of wave components, max- 
imum of 10 




16-25 


F 


Initial wave heading relative to 
craft heading (180° - head seas) 
deg . 


2, 3, 4.. 


. ’ 1-10 


F 


Wave length, ft. 




11-20 


F 


Wave amplitude, ft. 



Note: The program constructs irregular seas by adding to- 

gether a series of regular wave components with appropriate 
distribution of amplitude and frequency. To run a regular 
sea case, the number of wave components should be equal to 
1. If the number of waves is equal to zero (calm water) 
no additional data (cards 2, 3, 4...) should appear. 

Block 12 Initial Conditions 
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No option tags 



Card Columns Format 

1 1-5 I 

6-15 F 

16-25 F 

26-35 F 

36-45 F 



Entry 

01200 - control tag 
Initial speed, knots 
Initial pitch angle (positive- 
nose up) , deg . 

Initial draft at the center of 
gravity, inches 

Initial bubble gage pressure, psf 



Block 13 End of Case 
No option tags 

Card Columns Format Entry 
1 1-5 I 01300 - control tag 

Note: This card is used to signal the end of data for a 

given case. The program will now initialize all required 
variables and will start the calculation of the time history. 

If the program reads two consecutive 01300 cards with 
no data between, the job will be stopped rather than repeat 
the previous case . 



Block 14 End of Run 
No option tags 

Card Columns Format Entry 
1 1-5 I 01400 - control tag 

Note: This card is used to signal the end of the run. When 

encountered the program will print out the message: 

COMPLETED ALL RUNS 

and will then stop. 
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Block 16 Thrust and Rudder Input 



Option 


tag 1 = STB 


Thrust 


input 


Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01601 - Control tag and option 




6-15 


F10 


Value of constant thrust for the 
starboard screw if variable thrust 
is not desired 




16-25 


F10 


Number of data points in the 
thrust map. Left blank if 

constant thrust is used. 




26-35 


F10 


Side thrust coefficient 


2,3. . 


1-80 


8F10 


Data points for the independent 
variable, time in ascending order. 




. 1-80 


8F10 


Corresponding data points for 
the dependent variable. 



Option tag 2 = Port thrust input 
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Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01602 - Control tag and option 




6-15 


FI 0 


As above for the port screw. 




16-25 


F10 


11 11 11 li H || 




26-35 


F10 


M II 11 ll || || 


2,3. . 


1-80 


8F10 


11 II ll || m || 




1-80 


8F10 


11 H ll ll II || 


Option 


tag 3 = Rudder motion 


input 


Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01603 - Control tag and option 




6-15 


F10 


Value of rudder if constant rudder 








is desired 




16-25 


F10 


Number of data points of the 








rudder map. Left blank if 








constant rudder is desired. 


2,3. . 


1-80 


8F10 


Values of the independent variable 




• 




time in ascending order. 




’ 1-80 


8F10 


Corresponding values of rudder angle 








in degrees. Note: Positive rudder 








angle is right rudder. 


Block 


17 Bow and 


Stern Sea] 


Pressure Differences Input 


Option 


tag 1 = Bow 


seal 




Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01701 - Control tag and option 




6-15 


F 


Number of data points for the 



bow seal 
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Card 


Columns 


Format 


Entry 


2,3. . 


1-80 


8F10 


Time points in ascending order 


6,7. . 


1-80 


8F10 


Corresponding differences in 








bow seal pressures 



Block 17 shall be left blank if a constant pressure is desired. 
That constant shall be input in Block 5 and 6. However , . 
one seal may be held constant and the other varied. 

Option 2 = Stern seal 



Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01702 - Control tag and option 




6-15 


F 


Number of data points for the 








stern seal 


2,3. . 


1-80 


8F10 


As for option 1 , for the stern 








seal 


6,7. . 


1-80 


8F10 


As for option 1, for the stern 








seal 


Block 


18 Title 


card 




No option tags 






Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01800 - control tag 


2 


1-80 


20A4 


Alphanumeric data to be printed 








as title in input summary 


Block 


19 Fan Maps 




Option 


tag = 1 


Bow Seal 


Fans 
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Card Columns 



1 1-5 

6-15 
16-25 
26-35 
36-45 

2,3.. 1-80 



After 1-80 

last 

pressure 

data 

card 

Option tag = 2 
Card Columns 
1 1-5 ' 

• 6-15 
1.6-25 
26-35 
2,3.. 1-80 

After 1-80 

last 

pressure 

data 

card 



Format 

I 



F 



F 

F 

F 



8F10 



8F10 



Entry 

01901 - control tag and option 
Number of bow seal fans 
Speed of bow seal fans , rpm 
Number of data points (maximum of 25) 
If this is a 1.0 then the following 
data points are read. If left blank 
no further read takes place. 

Tabular value of fan pressure dif- 
ferential, psf. 

Corresponding volumetric flow rate, 
cu . f t . / sec . 



Cushion Fans 
Format 
I 
F 
F 
F 

8F10 

8F10 



Entry 

01902 - control tag and option 
Number of cushion fans 
Speed of cushion fans, rpm 
As option 1 above 
Tabular values of fan pressure 
differential, psf 

Corresponding volumetric flow rate. 



Option tag = 3 



Stern Seal Fans 



Card 


Columns 


Format 


Entry 


1 


1-5 


I 


01903 - control tag and option 




6-15 


F 


Number of stern seal fans 




16-25 


F 


Speed of stern seal fans, rpm 




26-35 


F 


Number of data points (maximum of 25) 




36-45 


F 


As option 1 above 


2,3.. 


1-80 


8F10 


Tabular values of fan pressure 
differential, psf 


After 

last 

pressure 

data 

card 


1-80 


8F10 


Corresponding volumetric flow rates 
cu . f t . /sec . 



Note: The tabular pressure data for all the fans must be 

stored in order of increasing pressure. (See discussion in 
Function FG1) . The data is for a single fan. 

Structure of Input Deck - The input deck for a given case or 
series of cases may be considered as being divided into two 
parts, "basic" data and specific data. Basic data is defined 
as those inputs which are required to adequately represent 
the craft regardless of the type of case being considered. 
Data such as mass properties, craft geometry, fan maps, etc. 
fall into this category. Specific data is defined as those 
inputs which are needed to define a specific problem such as 
runs in waves vs. calm water, rudder motion vs. no rudder 
motion, etc. 

The input cards which constitute the basic data are 
summarized below by control and option tag along with a 
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brief description of the input. An asterisk indicates that 
more than one input card is read in under the given control 
and option tag. 
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Required Input Summary 


(Listed in Numerical Order) 


Control Taq and Option 


Description 


00101 


Timing information 


00103 


Number of integrators and error 




tolerances 


00201 or 
00202* 


Mass and inertias or mass dstribution 


00300 


Craft geometric data 


00401 


Sidewall basic data 


00500 


Stern seal basic data 


00600 


Bow seal basic data 


00701 


Plenum geometry 


00702 


Froude number at hump speed 


00801 


Propulsion system basic data 


00901 


Rudder geometry 


01000 


Aerodynamic reference lengths 


01200 


Initial conditions 


01300 


End of data for one case 


01400 


End of all cases 


.01600 


Thrust map input 


01700 


Rudder motions input 


01901* 


Bow seal fan maps 


01902* 


Cushion fan maps 


01903* 


Stern seal fan maps 



*More than one input card under this control tag and 
option. 
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Optional Input Summary (Listed in Numerical Order) 



Control Tag 

and Option Description 



Program Action if Input 
Card Is Not Used 



00102 * 

00104* 

00105 



01101 * 

01800 



Subroutine No printout from individual 

print switches subroutine 

Tabular summary No tabular summary 

print switches 



Option swithces 
for lateral plane 
motions only, 
trim calculations 
and constant speed 
runs 



Six Degrees of freedom, fixed 
thrust (unless engine out is 
requested, surge velocity 
allowed to vary 



Wave input Calm water 

Title card No title printed 



Multiple Case Runs - In running multiple cases it is important 
that the user remember that any input data that is not read in 
again will be repeated from the previous case. Consider the 
case where a run is to be made consisting of a run in waves 
followed by a turn in calm water at a different speed. Assume 
both intermediate and tabular ouput is desired for both cases. 
In terms of control tags the deck might look as follows: 



00101 Required Input 

00103* 

00202 * 

00401 

00500 

00600 

*More than one input card under this control tag and option. 
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00701 



00801 




00901 




01000 




01200 




01600* 




01700 




01901* 




01902* 




01903* 




00102* 


Intermediate print 


00104* 


Tabular print 


01102* 


Wave input 


01300 


End of case 1 


01200 


New initial conditions 


01603 


Rudder motion 


01102* 


Wave input 


01300 


End of case 2 


01400 


End of run 



The purpose of the wave input card in case 2 is to read in 
the number of waves as zero (calm water) . Without this card 
case 2 would become a turn in waves. Similarly, if a run 
were to be made with case 1 being an engine out run and case 
2 being a turn with both engines operating, case 2 would 
have to reset the engine out time to a time larger than the 
run finish time by reading a large time on an 00802 card 

*More than one input card under this control tag. 
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(if case 2 had a finish time less than the engine out time 
from case 1 it would not be necessary to read this card) . 

Output - The output from INCON is printed after an end of 
case card (01300) is read. It is a combination of input 
data and calculated results. A sample case is given in the 
input at the end of the report. 
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(This page intentionally blank) 
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Subroutine INTGRL 



Discussion - Subroutine INTGRL is used to integrate a system 
of ifrst order ordinary differential equations. It uses a 
variable time step technique based on the Runge-Kutta-Merson 
algorithm. A fixed time step R-K also available - See Block 
1 option tag 3 for a discussion on its use. The program will 
stop the calculation if the time step becomes smaller than 
10 ® sec. 

Usage - CALL INTGRL (TIME) 

Descriptions of Parameters 

TIME = current value of the independent variable 

Input - The input to subroutine INTGRL is controlled by the 
00103 control tag as discussed in block 1 of subroutine INCON . 

Output - If the integration step size is halved, INTGRL will 
print the current value of time, the new value of the time 
step, the integrator number, the calculated error and the 
input maximum error as shown in the sample output at the end 
of the report. This output is useful as it tells the user 
which integrator is causing the step size to be reduced and 
may indicate an excessively tight input error tolerance 
(see discussion of 00103 card in discussion of subroutine 
INCON) . 
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Subroutine PROP 



Discussion - Subroutine PROP is used to calculate forces and 
moments on the craft due to the propulsion system. The sub- 
routine has the capability of uneven thrusting and receives 
this information by using the Function Program FGl to linearly 
interpolate a map input in block 16. 

Usage - CALL PROP 

Input - The input to subroutine PROP is controlled by block 
8 and block 16 of subroutine INCON and is discussed there. 

Output - The output from PROP is controlled by print option 
switch IPROP. If this switch is 1, the program will print 
the calculated forces and moments as shown in the sample 
output at the end of the report. 
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Subroutine RHS 



Discussion - Subroutine RHS is the subprogram containing the 

FORTRAN expressions for the right hand side of the system 

\ 

of first order differential equations. It contains the logic 
to calculate the bubble volume, area and pressure as a function 
of time as well as the logic which sums the individual forces, 
moments and leakages from the various craft components to get 
craft totals. RHS also contains the statements used for the 
writing of the scratch file for the tabular output. 

Usage - CALL RHS (VALUE) 

Description of Parameters 

VALUE = The array containing the values of the right 
hand side of the differential equations as 
calculated in RHS, i.e., for a system of 
equations of this form yy=fy(yp»t) the VALUE 
array is used to store the values of f^(yy^t) 

Input - None 

Output - The output from RHS is controlled by print option 
switch IHRS . When this switch is 1, the subroutine will 
print the output shown in the sample output in the back of 
the report. 
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Subroutine RUDDER 



Discussion - Subroutine RUDDER contains the logic for cal- 
culating rudder forces and moments as well as rudder motions. 
Rudder motions are calculated by employing the Function Pro- 
gram FGl to linearly interpolate a map input in block 16 of 
INCON. 

Usage - CALL RUDDER 

Input - The input to RUDDER is controlled by block 9 and 16 
in Subroutine INCON and is discussed there. 

Output - The output from rudder is controlled by print option 
switch IRUD. If this switch is 1, the program will print out 
the total forces and moments due to the rudder as shown in 
the sample output in the back of the report. 
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Subroutine SIDEWL 



Discussion - Subroutine SIDEWL calcualtes the forces and 
moments acting on the craft due to the sidewalls as well 
as the leakage flow rates associated with any gaps which 
open under the sidewalls due either to craft motion or 
waves . 

Usage - CALL SIDEWL 

Input - The input to subroutine SIDEWL is controlled by 
blocks 3 and 4 of subroutine INCON and is discussed there. 

Output, The output from SIDEWL is controlled by the print 
option switch ISIDWL. When this switch is equal to 1, the 
program will print the sidewall gaps, immersion depths and 
total forces and moments due to the sidewall-appendage 
combination as shown in the output listing at the back of 
the report . 
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Subroutine STNSL 



Discussion - Subroutine STNSL calculates the forces and 
moments acting on the craft due to the stern seal as well 
as the leakage flow rates arising from any gaps which open 
under the seal. The craft equilibrium leakage area is 
assumed to be located in the stern seal. 

Usage - CALL STNSL 

Input - The input data for STNSL is controlled by blocks 3 , 

5, and 17 of subroutine INCON and is discussed there. 

Output - The output from STNSL is controlled by print option 
switch ISTNSL. When this switch is equal to 1, the program 
will print the seal wetted lengths, and gaps as a function of 
station as well as the craft forces and moments due to the 
stern seal. An example is contained in the sample output at 
the back of the report. 
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Functions T1 and T2 



Discussion - Functions T1 and T2 are two small Function sub- 
programs used to calculate various trignometric relations used 
by subroutine WAVES . 
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Subroutine WAVES 



Discuss sign - Subroutine WAVES calculates the wave forces 
and moments acting on the craft. It also generates the 
wave heights at the various stations around the seals and 

sidewalls, as well as the bubble volume lost due to wave ele- 
vation . 

The sea state for irregular seas is computed in WAVES 
by adding together a series of regular waves with an ap- 
propriate distribution of amplitude and frequency. 

Usage - CALL WAVES 

Input - The input to subroutine WAVES is controlled by block 
11 of subroutine INCON and is discussed there. 

Output - The output from WAVES is controlled by print option 
switch IWAVES . When this switch is set equal to 1, WAVES 
will print the wave elevation around the seals and sidewalls, 
the wave dlevation at the center of gravity and the volume 
reduction of the bubble plenum due to waves as well as the 
total forces and moments acting on the craft due to waves. 
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USERS MANUAL FOR THE PLOTTING 
PACKAGE 

(INCLUDES TIIE MULTICURVE OPTION) 

DESCRIPTION: 

The plotting package of the COLFIL subroutine is a 
flexible routine that allows the user to: 

1) Plot on either the CALCOMP or the printer 

2) Plot up to 3 curves per graph on the CALCOMP 

3) Intermix multicurve CALCOMP and single curve 
CALCOMP plots 

4) Revise the variables plotted on successive runs 

5) Choose from up to 26 different variables 

6) Print up to 16 variables per run 

7) Revise the variables printed on successive runs 
The method of employing this routine must be adhered 

to or errors will result. JCL for calling this routine 
is found at the end of this manual . 

MEANING OF THE CARD 
1. Block 20 

Card Columns Format Entry 

1 1-3 13 020 

2 This card must be used if card 1 appears. See 
sample input 

1-10 II NCURV is the number of this 

curve on the corresponding 
graph, ie ; Col. 1 pertains to 
graph 1, Col. 2 - graph 2, etc. 
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Card Columns Format Entry 

A "0" must appear in these 
columns if only one curve per 
graph is to be plotted. "1" must 
appear if this is the first of 
more than one curve per graph, 
a "2" if this is the intermediate 
curve and a "3" if this is the 
last curve on this graph. 

2. Block 21 



Card Columns Format 

1 1-3 13 

2 1-16 812 

3 1-16 . 812 



3. Block 22 

Card Columns Format 

1 1-3 13 



Entry 

021 - Causes the summary variables 
to be read 

Punch in the numbers corresponding 
to the variable names desired to 
be printed in summary one. 

The same as card one but, for 
Summary Two, these cards need 
not be repeated for successive 
runs unless changes are desired. 

Entry 

022 - Read in the number of graphs 
desired, whether or not Calcomp is 
desired. Read the X and Y axes for 
the graphs and read the first line 
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Card 



2 



3 



Columns 



1-2 



3-4 



1-50 



1-48 



Format 



12 



12 



2512 



6A8 



Entry 

of the graph title. 

The number of graphs desired right 
justified, up to and including 10. 

A zero if print plot is desired or 
a 1, right justified, if CALCOMP 
is desired. 

The numerical values of the X and 
Y axes alternating in the following 
manner. The X axis of the first 
graph then the Y axis of the first 
graph and continuing for all 
desired graphs. A table relating 
axis number and name is at the 
end of this paper. 

The Alpha-numeric title that will 
appear on the first line of the 
CALCOMP plot. If CALCOMP is not 
desired, replace this card with a 
blank card. Your name should appear 
somewhere on this line for computer 
identification purposes. 
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SAMPLE INPUT DECK 

02Q 

11 This indicates that this is the first of 

multiple curves on graphs 1 and 2 . 



021 

0115 —a Summary one will print the first two variables, 

0116 ■ di, 15, and summary two the following two, 

01, 16. 



022 

0201 Two graphs are called for, 

0 1150116 , Graph 1 plots 01, 15, Time 

1 Graph 2 plots 01, 16, Time 



and a CALCOMP output 
and X displacement 
and Y displacement 



SES output for LEO 



Standard input deck 

013 



0201 
22 / 
021 \ 
1516' 
0115 



— As Above . 
runs . 



These can be changed on subsequent 



Changes desired 

013 

020 

33 

021 As Above. 

1325 
0116 



Changes desired 

013 

014 
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RELATION OF VARIABLE NAME AND NUMBER 



NUMBER 


NAME 


01 


Time 


02 


Wave Height 


03 


Z Displacement 


04 


Pitch Angle 


05 


Plenum Pressure 


06 


Bow Acceleration 


07 


C. G. Acceleration 


08 


Fan Power 


09 


Roll Angle 


10 


Yaw Angle 


11 


Lateral Acceleration 


12 


Speed through the water 


13 


Turn Radius 


14 


Turn Rate 


15 


X Displacement 


16 


Y Displacement 


17 


Air Flow in 


18 


Air Flow out 


19 


Net Force X Direction 


20 


Wave Force X Direction 


21 


Net Force Y Direction 


22 


Net Force Z Direction 
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NUMBER 



NAME 



23 


Net Torque X 


Axis 


24 


Net Torque Y 


Axis 


25 


Net Torque Z 


Axis 


26 


Rudder Angle 
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APPENDIX A-A 
CHANGING THE MODEL 
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Changing the Model 
Appendix A-A 

Introduction 

As previously stated, the program has been written 
in a modular form to facilitate changes to the program 
which reflect major design changes to the craft, i.e., 
a different type of seal, addition of controls, etc. Since 
it is obviously impossible to program all possible con- 
figurations in one program, these types of changes will 
require either replacing a whole module (sub-program) or 
altering the Fortran statements within a given sub-program. 

In order to achieve this, a more detailed knowledge of the 
structure of the program and the method of data trans- 
ference between subprograms is required. This appendix 
will attempt .to give a brief explanation of the program 
structure and will then discuss the details of how one 
would make certain anticipated changes . 

General Structure 

The primary variables (displacements, velocitites, pres- 
sure, etc.) are contained in an array of length 40 words, 
named VAL; which is transmitted through the various sub- 
programs using a labelled common block called VARBLE . Within 
each subprogram the various elements of VAL are equivalenced 
to local variables, which are given mnemonic names which 
closely resemble the standard mathematical symbols used for 

the variable (e.g. u, v, w, p become U,V, W, PB in the 

B 
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Fortran listing) . The second through the sixteenth elements 
of VAL are reserved for the output of integrators and must 
be used in numerical order, i.e., when adding a new equation 
to be integrated, the output of the integrator must be equi- 
valenced to the first unused location in VAL (VAL (12) for 
the current program) . 

The force and moment components acting on the craft are 
calculated for each element in the appropriate subprogram 
and are transmitted through various labelled cannon blocks. 
These force and moment components are then summed in sub- 
routine RHS to get the appropriate acceleration and velocity 
components. The accelerations and velocities are trans- 
mitted back to subroutine INTGRL in an array called VALUE. 
INTGRL then integrates these first derivatives to get the 
new velocities and displacements. 

Anticipated Changes 

a. Replacing Analytic Coefficients or Curve Fits With 

Tabular Experimental Data 

As experimental data is obtained it may be desired to 
replace coefficients based on analytic calculations or 
estimations from simple shapes by the experimental data. 

This data will often be tabular in nature. For functions 
of a single variable this is quite simple to do as the 
necessary table look-up routine already exists within the 
program (function FG1) . The user would be required to 
make the following changes. 
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1 . 



Statements to read in the tabular data must 
be added to subroutine INCON . 



2. A common block must be provided to transfer 



the data to the subroutine where it is to be 



used . 



3. A statement must be added to call FG1 to do 



the table look-up. 



A good example of the above procedure is afforded the 
user by the fan maps as currently programmed. This logic 
uses tabular data and the user should refer to input block 
19 of subroutine INCON and subroutine FAN to see an example 
of the three steps outlined above, 
b. Adding a Stern Seal Equation of Motion 

The current program has a stern seal model which is 
based on the assumption that the seal is completely de- 
formable (no stiffness) and as such there is no dynamic 
equation' of motion. For seal designs where there is an ap- 
preciable degree of stiffness a seal equation of motion 
should be added. 

Assume that the seal equation of motion to be added is 
a standard spring-mass-damper system and therefore may be 
represented by the following system of first order dif- 
ferential equations. 



G* = (-de c - ke + M )/m 
S s s 



( 1 ) 




( 2 ) 



171 



The following program changes would have to be made in 
order to incorporate Equations (1) - (2) into the current 
program. 

1. Logic must be added to block 5 of subroutine 
INCON to read in and store the constants d, 
k, and m as well as the initial conditions 
for 9 and 0 e . Care should be taken in the 
case of the latter to save these initial 
conditions using Fortran variables with names 
other than the names which will contain the 
current (running) values of the seal displa- 
cement and velocity. This is due to the fact 
that when running multiple cases the program 

will use the currently stored values of 9 and 

s 

9 S from the old case as initial conditions for 
the new case unless the initial conditions are 
reset. An example of this can be seen on page 
205 of the program listing after statement num- 
ber 1200 where the craft initial conditions 
are read in. Saving the initial conditions 
under a different name also has the advantage 
of allowing the input to be in degrees rather 
than in radians as required by the program, as 
the conversion can be done in the arithmetic 
replacement statements added in step 2 below. 
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2. In block 13 of INCON, (after statement 1300 
but before statement 1400) arithmetic re- 
placement statements should be added to con- 
vert the initial conditions to radians and 
radians/sec. and then store these values in 

the running variables for 0 and 0 , (see page 

s s 

205 of the program listing after statement 1302 
for an example) . Note that in the program 
the variable RAD is the number of degrees in 
a radian and may be used for the conversion 
of units (see example cited above.) 

3. The variable names for the seal, velocity and 
displacement should be equivalenced to VAL (12) 
and VAL (13) respectively. This automatically 
makes these variables the output of integrators. 
This equivalence statement should appear in sub- 
routines INCON, RHS and STNSL . 

4. The number of differential equations on the in- 
put card after the 00103 card should be increased 
from 10 to 12, indicating that two additional 
differential equations have been added. In 
addition, error tolerances for these integrators 
must be given (see page 154 of this report for 

a discussion of this point.) 

5. A common block to transmit d, k, and m must be 
added to subroutines INCON and RHS . 
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6. In subroutine RHS the Fortran equivalent of 
the right hand side of Equation (1) must be 
set equal to VALUE (11) while 6 S must be set 
equal to VALUE (12) . (See examples in listing 
of subroutine RHS.) 

7. In subroutine STNSL the expressions for M , the 

s 

total applied moment on the stern seal about its 
hinge point must be added. These are placed in 
STNSL as this subroutine has the logic for cal- 
culating wave elevations relative to the seal, 
etc. In addition the equations for the stern seal 
leakage gap is STNSL must be modified to take into 
account the seal motion. 

8 . A common block must be added to subroutines STNSL 

and RHS to transmit the value of the moment M to 

s 

RHS . 

c. Adding New Forces or Moments Acting on a Craft 

If the modification of the program entails the addition 
of any new force or moment component acting on the craft, care 
must be taken to add this component to the bending moment cal- 
culations as well as the craft motion calculations. If this is 
not done the equilibrium check done by subroutine SAM prior to 
calculating the bending moments will not balance and the cal- 
culation will be terminated. 

The following changes would have to be made: 
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1 . 



In subroutine RHS the new force component should 
be added to the WRITE statement preceding state- 
ment number 111 on page 218 of the program listing 
as well as into the calculation of GF in RHS. If 
the force component has a fixed moment arm, a 
common block should be added to transmit its value 
to subroutine SAM. If the moment arm is time va- 
rying the moment arm should also be added to the 
WRITE statement. 

2. In subroutine SAM the force and arm (if necessary) 
should be added to the following READ or WRITE 
statements; page 220 , statement numbers 42 and 
44; page 220 , after statement number 6. 

The force and moment components must be added 
to the equilibrium check on the top of page 
The force and moment components must be added 
to the calculation of the longitudinal and 
lateral moment and shear calculations. In 
these calculations the variables XMI(IX) and 
YMI(IY) are the distance to the current plane 
of interest. The user should study the current 
program listing (pages 183 - 233 to ascertain the 
procedure used to add these force and moment com- 
ponents . 



175 



APPENDIX A-B 

ABNORMAL TERMINATION MESSAGES 
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APPENDIX A-C 

SIDEWALL INTEGRAL TABLE PROGRAM 
DISCUSSION 
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Sidewall Integral Table Program 



Appendix A-C 

The sidewall Integral Table Program evaluates tables of 
integrals of sectional hydrostatic and hydrodynamic pro- 
perties over the whole and partial lengths of a SES side- 
wall. These tables are calculated as a function of sidewall 
draft, trim and length of a regular wave. The tables are 
then subsequently accessed by the main SES Motions and Load 
Program and when combined with knowledge of both the craft 
and waves instantaneous position and motion produces the 
hydrostatic and hydrodynamic forces and moments acting on 
the sidewalls. 

The primary input to this program is a table of water- 
lines and corresponding beams for each chosen station with 
arbitrary spacing in both waterlines and stations. The pro- 
gram is thus applicable to a general sidewall form. In 
addition to this input, the user must specify the range of 
drafts, trim (pitch) angles, and wavelengths (encompassing 
all components of a specified sea state) that the craft is 
expected to range thru in the subsequent simulations per- 
formed by SES Motions and Loads Program. The wave-lengths 
are specified via the wavelengths parameter (2TT/?0 • The 
first value of this parameter must be zero which corresponds 
to an infinite wavelength and is used to SIDEWL subroutine 
of the SES Motions and Loads Program to determine hydro- 
static and hydrodynamic forces due to body motion alone. 
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The tables at the remaining values are used by the WAVES 
subroutine for wave induced forces. Since the center of 
gravity is not necessarily known at the time this program 
is executed , the user must also specify a longitudinal 
reference coordinate measured from the transom about which 
the integrals are to be evaluated (a value near mid-ships 
is recommended) . The range of drafts and trims at the longi 
tudmal reference point must be selected in such a manner 



so that the top 


of the plenum does not sink below the water 


surface . 






The description and 


order of the cards of the input 


data deck is as 


follows : 




1st Card 






Columns 


Format 


Entry 


1-5 


I 


No. of wave pAfameters <20 


6-10 


F 


increment in wave parameters (ft. 


11-15 


F 


first value of wave parameter (ft 


16-20 


I 


no. of drafts < 5 


21-25 


F 


increment in draft (in.) 


26-30 


F 


first value of draft (in.) 


31-35 


I 


no . of trims < 7 


36-40 


F 


increment in trim (deg.) 


41-45 


F 


first value of trim (deg.) 


46-50 


I 


no. of stern beams <36 


51-55 


F 


increment of stern beam (ft.) 


56-60 


F 


first value of stern beam (ft.) 
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2nd Card 



Columns 


Format 


Entry 


1-5 


I 


no. of stations = NSTA <95 


6-10 


I 


max. no. of drafts per station <8 


3rd Card 






Columns 


Format 


Entry 


1-5 


F 


distance from keel to top of plenum (ft.) 


4th Card 






Columns 


Format 


Entry 


1-10 


F 


distance of reference position 






from transom (ft.) 


5th Card 






Columns 


Format 


Entry 


1-5 


I 


no. of integral limits in addition 






to the integral over the whole length 2,10 


6-75 


10F5 


integral limits (distance from transom) 






(ft.) 


Next NSTA Cards 


(one card 


per station) 


Columns 


Format 


Entry 


1-5 


F5 


distance of station from bow (in.) 


6-10 


F5 


lowest waterline beam (ft.) 


11-15 


F5 


corresponding height of beam (ft.) 


16-20 


F5 


next waterline beam (ft.) 


21-25 


F5 


corresponding height of beam (ft.) 
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Columns 



Format 



Entry 



66-70 


F5 


next waterline beam 


(ft.) 


71-75 


F5 


corresponding height 


of beam (ft.) 



This program transmits the calculated tables to the SES 
Motions and Loads Program by writing the tables on units which 
are later read by the main program. Unit number 10 contains 
the tables obtained by integrating over the entire length and 
the succeeding units contain the tables obtained by integrating 
over partial lengths . 
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APPENDIX B 

FORTRAN PROGRAM LISTING 
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